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On the formation of shocks of electromagnetic plane waves in non-linear 

crystals 


Demetrios Christodoulou and Daniel Raoul Perez 


Abstract. An influential result of F. John states that no genuinely non-linear 
strictly hyperbolic quasi—linear first order system of partial differential equations 
in two variables has a global C 2 — solution for small enough initial data. Inspired by 
recent work of D. Christodoulou, we revisit John’s original proof and extract a more 
precise description of the behaviour of solutions at the time of shock. We show that 
John’s singular first order quantity, when expressed in characteristic coordinates, 
remains bounded until the final time, which is then characterised by an inverse 
density of characteristics tending to zero in one point. Moreover, we study the 
derivatives of second order, showing again their boundedness when expressed in 
appropriate coordinates. We also recover John’s upper bound for the time of shock 
formation and complement it with a lower bound. Finally, we apply these results 
to electromagnetic plane waves in a crystal with no magnetic properties and cubic 
electric non-linearity in the energy density, assuming no dispersion. 


1. Introduction 

Electromagnetic plane waves of a single polarisation can form shocks in non linear 
media ILL84I §111]. If both polarisations are present, one may appeal to I.Ioh74l in the 
anisotropic case. A solution to the general three-dimensional problem of electromagnetic 
waves in non-linear media, even in the homogeneous and isotropic setting, is still out of 
reach. A model case, however, was recently dealt with by Miao and Yu IMY14| . 

In this paper, we focus on electromagnetic plane waves in non-linear crystals. These 
are solutions of certain quasi-linear hyperbolic systems of first order partial differential 
equations, which fall in the framework of Fritz John’s treatment |.Toh74l . However, the 
description in I )■ 'Ii7 II is somewhat incomplete, since only the non-existence of smooth 
global solutions is shown. The purpose of the present work is to overcome this limitation. 
Its key feature is a more precise description of the behaviour of smooth solutions at 
their final time of existence. It turns out that their singular behaviour, as mentioned 
by John, is a consequence of the singular relationship between the physical space-time 
coordinates and coordinates adapted to the characteristics. Indeed, all relevant quantities 
remain bounded when expressed in appropriate coordinates adapted to the characteristics, 
and a shock is characterised by an inverse density of characteristics decreasing to zero 
somewhere. This is then similar to the situation in |Chr07l . where the formation of shocks 
in relativistic compressible fluids is considered (but see also the more recent and self- 
contained exposition of the non-relativistic case ICM14I 1. Those shocks are characterised 
by the foliation into outgoing characteristic null-hypersurfaces becoming infinitely dense 
somewhere. At the same time, the solution can be smoothly extended to the time of shock 
when expressed in acoustical coordinates. This description served as an inspiration for 
the work at hand. 

Our results are not limited to the four-by-four system obtained when considering both 
polarisations of electromagnetic plane waves in non-linear crystals, but they apply to the 
general framework set up in lJoh74l . i.e., genuinely non-linear strictly hyperbolic quasi- 
linear first order systems of partial differential equations in one space variable and one time 
variable. One calls a system strictly hyperbolic when the characteristic speeds are distinct, 
whereas such a system is said to be genuinely non- linear, if certain directional derivatives 
of the characteristic speeds do not vanish. As in IJoh74l , a key ingredient is the assumption 
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of smooth, compactly supported and small non-zero initial data, leading to the long-time 
behaviour of solutions as essentially simple waves, i.e., maps with one-dimensional range. 
This is expressed by the fact that the problem reduces to a system of non-linear ordinary 
differential equations along each characteristic, obtained by an appropriate decomposition 
of the spatial derivatives of the solution vector. Calling characteristic strip the collection 
of characteristics of a certain speed emanating from the support of the initial data, John’s 
result states that solutions have to remain (7 1 -bounded outside those strips. On the other 
hand, John shows a blow-up in finite time of a component of the first spatial derivative of a 
solution inside at least one of those strips (although solutions themselves remain bounded). 
His argument can be modified slightly to get a deeper insight into the failure of regularity. 
By expressing the solution in characteristic coordinates rather than the Euclidean ones, 
the same component of the first spatial derivative remains bounded in every strip. This 
is due to the fact that, in those coordinates, the spatial derivative corresponds to the 
partial derivative with respect to Euclidean space multiplied by the inverse density of 
characteristics of the strip at hand. This inverse density acts as a regularising factor. The 
final time is then characterised by this quantity going to zero somewhere, which, therefore, 
implies John’s result. 

The first part of our work can, in fact, be extracted directly from John’s proof. John 
introduces the inverse density of characteristics, as well as the regularised component of the 
first order spatial derivative of the solution. However, he uses both those quantities only 
in an integral sense, thereby hiding their actual behaviour at the final time. This is why 
we choose to re-expose completely his proof. We also use slightly different assumptions on 
the initial data. While he requires continuous second derivative, we only need Lipschitz- 
continuity of the first derivative. Even if this is not an important modification, we feel 
that it is more natural, as the initial system is only of first order. 

In a second step, we extend the above result by considering the derivatives of second 
order of solutions. It is not hard to see that, again, outside every characteristic strip, all 
quantities of interest remain bounded. Inside such a strip, however, the situation becomes 
more complicated, as the quantities relating to the strip at hand depend too strongly on 
quantities related to other strips. The problem is that we cannot bound certain quant¬ 
ities when integrating along the appropriate characteristics, while inside a characteristic 
strip of a different speed. This is overcome by expressing the relevant quantities in bi¬ 
characteristic coordinates, where time as a parameter for a characteristic is replaced by the 
spatial characteristic coordinate of a transversally intersecting characteristic of a different 
speed. 

We end this paper with an application of the theory to electromagnetic plane waves 
inside crystals with no magnetic properties and cubic electric non-linearity in the energy 
density, neglecting effects of dispersion. One should point out here, though, that physical 
experiments would always exhibit dispersion which will have the tendency to counteract 
shock formation. Accepting that a description of the full three-dimensional problem in 
crystals is still far away, solving the present problem seems like a good complement to a 
potential three-dimensional theory for isotropic media. Indeed, one would probably have 
to focus the beam in the latter case in order to counteract dispersion, as the non-linearity 
is only of fourth order in the Lagrangian. 

The present work is organised as follows. In Section [5] we introduce the necessary 
notation in order to state our result. We explain characteristics and the coordinates 
related to them. We also expose some first order variations and formulate the technical 
assumption for hyperbolic systems to be genuinely non- linear. Section [3] is devoted to 
a re-exposition of John’s proof. We begin by obtaining the evolution equations along 
characteristics of first order quantities. After presenting some a-priori assumptions on the 
initial data, we recover John’s bounds on those quantities using his method. We end the 
section by analysing more in depth the behaviour near the final time of the inverse density 
of characteristics. In doing so, we establish upper and lower bounds on the final time (but 
see also IHor97l for existence, uniqueness and lifespan of solutions). Moreover, we show a 
lower bound on the inverse density of characteristics in certain regions, including parts of 
its relative strip. In Section|4j we prove boundedness of the second order quantities. Using 
characteristic coordinates for a particular characteristic strip, the respective quantities also 
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depend on second order quantities related to other characteristic speeds. We deal with 
the latter by expressing them in bi-characteristic coordinates. It turns out that they form 
a complete regular system that can thus be bounded. The bounds on the second order 
quantities related to the strip at hand are then obtained by a Gronwall-type argument 
for systems. We end this paper with Section [5j in which we consider electromagnetic 
plane waves inside crystals with no magnetic properties and cubic electric non-linearity 
in the energy density, assuming absence of dispersion. After quickly recalling the relevant 
equations, we start by exposing under what conditions the theory can actually be applied. 
We then point out that the simplest two-polarisation case actually decouples and can be 
dealt with using Riemann invariants. In contrast, as we subsequently show, no Riemann 
invariants exist for generic cubic non-linear energy densities. As an actual experimental 
setup would require to deal with a material boundary, we end this section with giving a 
sketch of that situation. 


2. Preliminaries 


2.1. The Cauchy problem. Let TV £ N \ {0} and U C M N be open with 0 £ U. Let 
K £ N\{0} and a £ ( 7 max { 2 > x 1 ([/; M. NxN ) be such that the eigenvalues Xi (i £ {1,..., TV}) 
of a fulfil Ajv < ... < Ai in U, i.e., we assume that the first order quasi-linear system 

\ <9 , . d 

(!) _ M + a(M )_ M = 0 

is strictly hyperbolic around the trivial solution u = 0. Take S > 0 small enough so that 
the open TV-ball T?^(0) of radius 2 S around 0 £ is completely contained in U, i.e., 

(2) B? s ( 0) C U. 

For each i £ {1,..., TV}, define 

(3) Xi = sup Xi(v), = inf A i(v), 

and 


(4) a = min(A fc — A;). 

k,l 

k<l 

We will assume 5 to be small enough so that 

(5) a > 0, 

which is possible, since a is continuous. This means that the system 0 is uniformly 
strictly hyperbolic in B£s( 0). 

For a general TV-dimensional vector space V endowed with a basis {T5i}ig{i,...,iv}, we de¬ 
note by {E* x } ie { the basis of the dual space V* of V which is dual to {£'i}i6{i,...,iv}, 
i.e., we have 

E* i Ej — 5 i j, Mi,j £ {1,..., TV}. 

We can then expand any vector V £ V and any covector r/ £ V* in the bases {} ie {! ... ,rv> 
and {F ; *‘}i6{i,...,iv}, respectively, as follows: 

V = J2(E* i V)E i , r, = ^2(r)Ei)E*\ 


We define the scalar product (■ , • )e on V by 

{V,W) E = J2(E* k V)(E* k W), V,WeV. 

k 

Then 

(Ei,Ei)B = Y, 5k ^ k i = M € {1,..., TV}, 

k 

so that {-ET}i6{i,...,jv} is an orthonormal basis of V with respect to (•, - )e- Let le denote 
the canonical isomorphism of V* onto V induced by (•, - )e- It is defined as follows. Given 
a covector r] £ V*, the element le (i?) £ V is the unique vector which fulfils 

(r E (v),V) E = V V, MV £ V. 
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It follows that ie maps E* 1 to Ei, i £ {1 Also, the isomorphism le induces a 

scalar product (•, -)e* on V*, defined by 

(v,Oe* = »7i£eV*, 

relative to which the basis is orthonormal: 

(E*\E* j )E+ = < t E (E* i ),L E (E* j )) E = (E i ,E j ) B = S ii , i,j £ 

Returning to the system ©. let eoi = ei(0) denote an eigenvector of o(0) correspond¬ 
ing to the eigenvalue Ai(0), i £ {1,..., N}. Since all the eigenvalues of a are distinct, 
{eoi}ig{i,...,jv} forms a basis of M. N . Denoting by {eo**}ig{i ! ... ! jv} its dual basis, we define 
the scalar product (•, ■ )o by 

(6) (v, w}o = ^(eo* fe «) (eo* k w), v,w£R n . 

k 

With to the isomorphism of R^ onto M. N induced by (•, • )o, the inner product (•, ■ )o* 
induced by to on R N is given by 

(t7,£)o* = (to(»7),to(C)) 0 J ri,^eM. N . 

We then obtain the induced norms || • ||o and || • ||o* on R lv and R^*, respectively, from 

(7) || V || 0 = \J(v,v)o, v £ R n , and ||r?||o* = \J{v,v)o*, V £ R N *. 

Now, for i £ {1, ...,N} and u £ B% s ( 0), let et = efiu) be the eigenvector of a(u) 
corresponding to the eigenvalue A i(u) and fulfilling the normalisation condition 

(8a) ||ei|| 0 = 1, * £ {1,..., N}. 

By the strict hyperbolicity of ©, the ei form a basis of R w . Let {e* I }ig/i i ... i jv} denote 
the basis dual to {ei}j 6 .n so that 

(8b) e* l e j =8 l j , i,j £ {1,, N}. 

Observe that a, and thus e* 1 , are determined uniquely up to sign by (I8al) . Moreover, 
viewing a as a linear map acting on the dual R^ of R w , we have 

(9) e* l a = \ie*\ i€{l,...,N}. 

Notice, as well, that with the matrix a, also a, e** and A i are ( 7 max f 2 . K I 

Remark. Defining the (N x N)-matrices G and Go by Gij = d ■ e-j and Goij = {ei, efio, 
i,j £ {1,..., N}, respectively, e* 1 , i £ {1,..., N} acts on vectors v £ R^ by 

e* l v = ^(G _1 ) IJ ej • v = '^{Go^ 1 ) 13 {e j ,v) 0 . 

3 3 

Here, denotes the Euclidean inner product in R N . The normalisation condition (18a.lt 
states that Goa = 1, i £ {1 This is different from I.Ioh74l . where John sets 

(G- l )« = l, i £{!,...,JV}. 


Remark 1. For a better geometric understanding, consider the following. Think of u as 
a mapping u : A4 —> M, where M = R x [0, T] is (a strip in) the ( x,t)-plane and A f is 
an N-dimensional manifold. Then, given {xo,to) £ M and writing yo = u(xo,to) £ M, 
the vectors -^u{xo,to) and -^u(xo,to) in T yo Af are tangent to the images under u of the 
coordinate lines {t = to} and {x = xo\, respectively, in the ( x,t)-plane. Moreover, a(yo) 
is a map of T V0 N into itself, and the equation JT} at (xo,to) gives a relation between the 
tangent vectors u(xo,to ) and f^u(xo,to): 



+ a(u{xo,t 0 )) ^^M(*o,to)^ 


= 0 . 


A smooth vectorfield v on M is an eigen-vectorfield of a, if v(y) is an eigenvector of a(y) 
at each y £ A f. So there is a smooth function A on Af such that 


av = \v. 


Note that the condition that the eigenvalues X of a be real and distinct implies that the 
manifold A f is parallelisable. 
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At any arbitrary y £ M, we can also consider a(y) as a linear map from T*J\f into 
itself by: 


c H* (a(y), VC e T*M, 


where 

{Ca(y))w = C(a(y)w), Mw € T y Af. 

A smooth 1-form on M is an eigen-l-form of a, if £(y) is an eigen-covector of a(y) at 
each y £ A f. So there is a smooth function X 1 on A f such that 

& = X'C 


Assuming all the eigenvalues X j, i £ {1, ...,AT} ; of a are distinct, and given a basis 
( e *}i6{ 1 ,...,Jv} °f eigen-vectorfields, the dual basis {e* t }i e r 1 _ N \ of 1-forms fulfils 

e *\ y ) ej ( y ) = S l j , Mi, j £ (1,..., N}, My£Al. 

It follows that e* z is an eigen-l-form of a with the same eigenvalue function Xi as e*. 

Now consider a non-zero “seed” function /o € C^’ 1 (R; 77^(0)) with support in the 
interval 7o = [—1,1]. Here, C^’ 1 is the space of compactly supported K times differentiable 
functions with Lipschitz continuous K th derivative. For ■& > 0, let 


( 10 ) 


/ = 0/o. 


In what follows, we want to study solutions u £ C K (R x [0, 7*); 73^(0)) of the Cauchy 
problem 


( 11 ) 



( x , t) + a[u{x, t)) -rp-u(x, t) = 0 


u(z, 0) = f(z), 


(x,t) eRx [0,7*), 

zel, 


on their maximal time-slab of existence R x [0,7*). Our focus will lie on their behaviour 
at the final time 7*. 


2.2. Characteristics (cf. f Joh740 . Consider a solution u £ C k (R x [0,T];73^(0)) of 
d) for some T > 0. For z £ R and i £ (1,..., N}, let Ci(z) denote the 7 th characteristic 
starting at z, i.e., the set of points 


( 12 ) 


Ci{z) = {(Xi{z,t),t) elx [0, T\ 



where Xi(z,t) is a solution of 


(13) 


^- f Xi(z,t) = Xi(u(Xi( Z ,t),t)y 
Xi(z, 0) = z. 


Note that, since A i,u £ C K , the theory of ordinary differential equations tells us that X t 
is unique and itself C K . Therefore, 


(14) Xi(zi,t) < Xi(z 2 ,t), for all 7 £ [0,T], whenever zi < Z 2 , 

so that 


(15) Xi (—1,7) < Xi{z,t) < Xi{l,t), for all z £ Iq = [—1, 1], 7 £ [0,T[. 

Notice that u is supported in the region delimited by 7 = 0, 7 = T, and the extreme 
characteristics Cn(— 1) and Ci(l), which are actually straight lines: 

(16) A'jv(—1,7) = —1 + Ajv(0)7, Xi (1,7) = 1 + Ai(0)7. 

We now show that our assumption (0 ensures that the sets Hi = Ci(I 0 ) = U ze / 0 c 4 2 )- 
which we shall call characteristic strips, are “separated” after a short time. Indeed, from 
CD, ©), we have, for all* € (0,..., N}, z £ Io, 7 £ [0, T], 

(17) z + XJ < Xi(z,t) < z + Xit. 

But then, using 0 , we have for k < Z, 

Xk(— 1, t) — Xi (1, t) > (—1 + — (l + Ajt) ^ —2 + (A fc — Xi)t > —2 + crt, 
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which, by ©, is strictly positive as soon as t exceeds 

(18) to = -. 

G 

As a result, 

(19) (7l fr p)7Ji)p){(a:,t) | te [*0,211=0, Vfc + l. 

Next, we observe that m implies that we can introduce a new set of coordinates on 
R x [0, T], 

2.3. Characteristic coordinates. Let i £ {1,..., IV} be given. Since, by (1111) . Xifzi. t) < 
Xi(z 2 ,t), whenever zi < Z 2 , and since Xi £ C K , we can introduce the following coordin¬ 
ates ( Zi,Si ) on R x [0,T]: Given ( x,t ) £ R x [0,T], there is a unique ( Zi,Si ) £ R x [0, T] 
such that 

(x,t) = ( Xi(zi,Si),Si ). 

Introducing the function (cf. [Jo m) 

( 20 ) pi(zi,Si) = -^-Xi(zi,Si) 

UZi 

on R x [0, T], and recalling from (1131) that 

\ i ( z i , s i ) = £- Xi ( Zi , s i ), 

we obtain 

( dx = pidzi + A idsi, 

1 dt = dsi. 

Note that pi represents the inverse density of the I th characteristics, whereas A i rep¬ 
resents their speed. It will turn out that pi acts as a regularising factor for u inside the 
characteristic strip IZi. The formation of “shock” will then be characterised by pi tending 
to zero. 


! _d_ _ d_ 

dzi ~ pl dx ’ 

_d_ _ d_ d_ 

dsi 1 dx + dt 


2.4. Bi—characteristic coordinates. Observe that, by our assumption ©, the inter¬ 
section between an I th and a j th characteristic (i ^ j ) has to be transversal. Fixing then 
i,j £ {1,..., N}, i ^ j, we can eliminate Si in favour of Zj as a coordinate, Zj acting as 
a parameter along the characteristic Ci(zi) (and likewise for Sj). The coordinates thus 
obtained on R x [0,T] are described as follows: Given ( x,t ) £ R x [0, T], there is a unique 
( Ui-iVo) ^ such that 

(M) = (Xi(yi,t'(yi,yj)),t'{yi,yj)^ = (^X^y^t'{yi^j)) ,t'(yi, y^, 
for some C K -function t' of yi and yj. We compute 



It follows that 

—t' = pi —t' = Pj 

dyi A j A, Oijj A, A j 

where the denominators are non-zero by ©■ Hence, using also ED, 


(22a) 


d _ pj\j d pi d _ pi d _ _d_ pi d 

dyi A j — A i dx Xj — Xi dt X j — X i dsj dzi Xj — Xi dsi ’ 

9 _ PjXi d pj d = pj d = _d_ pj d 

dyj Xi — Xj dx Xi — Xj dt Xi — Xj dsi dzj Xi — Xj dsj ’ 
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and 


(22b) 


dx = 

dt — dsi — ds j — 


PAj j.. | Pj Ai 


-dyi + 


Ai - A. 


■dyj , 


~^~T dyi + \ Pj \ dyj ’ 

i Aj Ai Ad 


dzi = dy t , 
dzj = dyj. 

Before stating our results, we introduce the technical assumption ultimately responsible 
for the formation of shocks. 

2.5. The first variation of a and genuine non-linearity (cf. |.Toh74n . So far, we 
only required our Cauchy problem (Hill to fulfil ([5]), i.e., we assumed the matrix function 
a of u to be strictly hyperbolic around the trivial solution in such a way that the different 
characteristic speeds are uniformly bounded away from each other. In order to ensure the 
formation of shock, we require, in addition, that the system ID be genuinely non-linear 
around the trivial solution, i.e., for each i, we assume that the change of Ai(0) in the 
direction ei(0) is non-zero. 

Before giving the precise definition at the end of this subsection, we introduce the 
following coefficients to describe the first variation of a: 

(23) 4z = t 4z( u ) =e* : ’(u)(D ei(u) a(u))e k (u), j, k, l € (1,..., N}, 

where 


(24) 


D, 


e z0 


i)<p(u) = dip(u)ei{u) = lim - ^ p(u + eez(«)) — ip(u)^j 


denotes the directional derivative of a (scalar-, vector- or matrix-valued) function ip of u 
in the direction ez. Note that, since {ei(w)}-. 1 N forms a basis of the tangent space 

T U M. N at u, the coefficients defined above are sufficient to completely describe the first 
variation of a. 

Now, since e k is an eigenvector of a, so that (a — A*, id)e*, = 0, we have 
(a - \ k id )(D ei e k ) = ~(D ei a - D e ,\ k id)e fe . 

Thus, using (f8bl) and (ED, we get 

(A m - A k )e* m (D El e k ) = e* m (a - X k id )(D e ,e k ) 

= —e* m (D ei a — D ei \ k id)e fe = -c kl + S rn k (D ei \ k ). 

Setting m = k, this yields 


(25) 

while for m ^ k we get, 


D ei Xk — Ckl , 


■km / 7 -v \ 

e (D ei e k ) = 


Xk X r 


~^kl • 


Thus, expanding D ei ek in the basis {ei,..., ejv}, we have 

Dei e fc — ^ ^ ( e (D ei e k )^em = (e (D ei e k )^e k + ^ ' t , 

— — Ak An 


'Ckl€"n 




But, since ||efc||o = 1, whence {{D ei e k ),e k \ = \ D ei (||e fc ||§) = 0, we obtain 
e* k (De l e k ) = — V v— ^——cZ(em,e k ) 0 , 

‘ ■* A l- - Am 


from which we infer 


(26) 


Dei Cfc — 


‘ ^ Xm Xk 

my£k 

From the duality condition (|8bll. we get 


c£z((e m 5 e k )oe k — Cmj • 


(■ De l e k3 )e k = -e* J {D ei e k ), 
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and compute 


{D ei e* 3 )e k = - ^2 3 c£)((e m ,e fc )o<4 ~4n)- 


m^k 


Thus, expanding D ei e* 3 in the basis {e* 1 ,..., e* N }, it follows that 

D ei e *° = '^2{{De l e* J )e k )e* k = ^ x \ 4 (( e m, e k ) 0 5 J k - 5 ] rn )e 


a^k 


= E J-rxZ cTl{ern ’ ei)oe * j ~ E 




k 

k^j 


(27) 


= E A. - A fc ( c ^ efc ’ e ^° e * J +c Li e * k )- 


k 

k^j 


We are now ready to give the definition of genuine non-linearity for the system Q: 
we simply require c^(0) = D ei \i(0) ^ 0 (cf. (I!?51) and (ESI) ). 

By an appropriate choice of sign for each ei (and hence for e* 1 ), we can without loss 
of generality assume that 

(28) 4(0) <0, V* € {1, ..., IV}. 

Choosing then <5 > 0 small enough, we can ensure by continuity that 

(29) 44X0, V« £ Bm(0), Vi € {1,..., N}. 

2.6. Statement of results. Before we can state our results, we need to introduce the 
following first and second order quantities related to a solution u £ C K (iBi x [0, t*); B^( 0)) 
of the Cauchy-problem m subject to the additional assumption (1281) on its maximal 
time-slab R x [0,t*) of existence. 

Fix i £ {1,..., N}. For k £ (1,..., N}, let 

(30) w k {x,t) = e* k (u(x,t))-^-u(x,t), 
and, for the given i, let 


(31) 


\zi,Si ) = e* l (u(zi,Si))-^-u(zi,Si) = pi(zi,Si)w\zi,Si). 


dzi 


Here, and in the following, we abuse notation and write 

u(zi, Si) for u(Xi(zi, Si), Si), etc... 

Further, define 


(32) 


d 


d 


Pi(zt,Si) = —pi(zi,Si) and v\zi,Si) = —v\zi,Si). 


dzi 


dzi 


Finally, fix j £ {1,..., N}, j ^ i, and define 
d 

(33) (l 4 (y t , yj) = — pj {yi, yj ) and 

oyj 

Theorem 2. Let N, K £ N \ {0}, let U C R N be an open neighbourhood of 0, and let 
a £ (7 max l 2 ’ K l ([/; R iVxJV ) be such that the system 


( ’X(i n,yj) = — v 3 (yi,yj). 
oyj 


CD 


d . , d 
_ u + a(M) _ M = 0 


is strictly hyperbolic and genuinely non-linear. Let 8 > 0 be sufficiently small, so that 
conditions 


m 

© 

and 

dHJ 


Bm(0) c u, 

a > 0 

44) <0, Vu £ Bm(0),« £ {!,•••, N}. 







On the formation of shocks of electromagnetic plane waves in non-linear crystals 


9 


hold. Let fo £ C^’ 1 (R; B^ (0)^J have support in Iq = [—1, 1]. Finally, let e £ (0, y-hj) be 
a parameter. 

Then there is a constant i9o > 0, depending solely on a, 6, e and the C K ’ 1 -norm of fo, 
such that the following statements hold for any solution u £ C K ( R x [0, £*); B^(0)^ of 


dJ 



{x,t) + a(u{x,t )) — u(x,t ) = 0 


u{z, 0) = /(z), 


(x, t) £ R x [0, t *), 
z £ R, 


with initial data f = flfo and given on its maximal time-slab of existence R x [0, t*), 
whenever 0 < i? < i?o- 

0) (cf-\ 2 oh74l ) For every i £ {1, ..., N}, we have: pi and u* are bounded everywhere. 
Moreover, pi(zi,t) > 0 for all Zi £ R, t < t *. Also, Wi is bounded outside the strip 

7 Zi. 

(ii) Let Wq = sup w k (z, 0). Then 

fee{ 

z£R 


min 

i 


_1_ 

(l + e) 3 |4(0)|W 0 + 


< t* < max 

i 


_l_ 

MW?' 


Moreover, there is an index i + £ (1,..., N} and an initial location z + £ R such 
that Mm.t/'t * pi+(z + ,t) = 0. 

(iii) If K > 2, we have for every i £ {1,..., N}: pi and v l remain bounded inside IZi. 
Moreover, for every j £ {1,..., N} \{i}, and remain bounded inside TZi. 
Also, -§^u> 1 is bounded outside all the strips, i.e., in (R x [0, £*)) \Ufeg{i jv} 


2.7. Some notation. In the following and similarly to IJoh74| . we shall write A = O(B) 
for two quantities A and B depending on the initial data / = dfo, whenever there exists 
a constant C, depending only on the matrix function v i-> a(v), on the size 5 > 0 of the 
IV-ball on which our system © is uniformly strictly hyperbolic, in the sense that © 
holds, as well as genuinely non-linear, in the sense that CS0> holds, and on the C K,1 -norm 
of the seed fo, such that 

\\A(f)\\ <C\\B(f)\\, 

for all / with i9 > 0 small enough. Here, the meaning of || • || depends on the quantity under 
consideration. It is the usual absolute value when applied to scalars. If the quantities are 
vectors, we use the Euclidean norm | ■ |, the norm || • ||o or the norm || • ||o* (cf. 0), in 
dependence of whether the quantity lives in B™( 0) C R w , in the tangent space T U R N or 
in the co-tangent space T^R^, respectively. 

Also, for the rest of the paper, let e £ (0, ) be fixed. We will use it as a parameter 

to tweak our estimates so as not to get too big deviations, the idea being that e is very 
close to zero (assuming e < 1/100 is enough, but we can obtain better approximations if 
we choose it even smaller). 


2.8. A reminder on the Gronwall lemma. For the sake of completeness, we recall the 
precise statement of the Gronwall lemma in the form in which we are going to apply it. 
We do this to avoid stronger regularity conditions or additional assumptions on the signs 
sometimes encountered in the literature. Note also that the lemma remains valid with 
reversed inequalities. 

Lemma 3. Let T > 0 and let rj be an absolutely continuous, real-valued function on 
[0,T). Also, let if and ip be two locally integrable, real-valued functions on [0,T). Assume 
that 

(34) ^ V(t ) < ip(t)v(t) + ip(t) 

is satisfied for almost every t £ [0, T). Then the following inequality holds for all t £ [0, T) 
rj(t) < Jo ^ r)dr ^(0) + J* e~ J ^ r)dr <p(s)ds 


( 35 ) 
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Proof. By the fundamental theorem of calculus, we have -d- J 3 ip(r)dr = ip{s) for almost 
every s £ [0,T), since if is locally integrable. Therefore, for almost every s £ [OjT), 

^(v(s)e~S > <* e~So*^%( s ). 

Integrating this from 0 to t £ [0, T), we obtain for every t £ [0, T) that 

v(t)e~ So r)dr - n{ 0) < f e~ So ' Kr)dr <p(s)ds, 

Jo 

from which ([35]) is immediate. □ 

We deduce the following result for a special case of two coupled integral inequalities. 

Lemma 4. Let T > 0 and let ?;o, */i be two non-negative, locally integrable functions on 
[0, T). Assume further that there exist non-negative constants Ao, A\, Bo, B\, as well as 
two real-valued, locally integrable functions po, Pi on [0, T), such that 


(36) 


Vi(t) < Ai / rn(s)ds + Bi / r}i-i(s)ds +i £ {0,1}, 


is satisfied for almost every t £ [0, T). Then there are constants cl, d, f J , gt (i,j £ {0, T}), 
depending only on Ao, A\, Bo, B\ and T, such that the following inequalities hold for every 
t £ [0, T) and every i £ {0,1} 

(37) rji(t) <^2ci (d ( e d( * _s) f e fJ( - 3 ~ r) pj(r)drds+ j e 9i ^~ s) pj(s)ds J + Pi(t). 

j=0 \ ■'O "'° / 

In particular, if po and pi are bounded on [0, T), then so are rjo and r]\. 

Proof. For t £ [0, T) and i £ {0,1}, let 

Vfit) = / rn(s)ds. 

Jo 

Then ip is absolutely continuous, hence almost everywhere differentiable, as well as locally 
integrable. Moreover, rji is non-negative, non-decreasing, and we have rji(0) = 0. From 
dm), we see that for almost every t £ [0, T), 

^rj-fit) < Afifjfit) + Bifji-i(t) + pfit). 

Then we obtain from Lemma [3] that, for every t £ [O,! 1 ) and i £ (0,1}, 


»}<(*) < 
Letting 


Ai(t-s) 


(Birn-i(s) + pi(s))ds < e AiT Bi j r)i-fis)ds + f e Ai(t 3 ^pi(s)ds. 

Jo Jo 

pfit) = f e Ai(f ~ s ' > pfis)ds, 

Jo 


Bi = e AiT Bi 


we can write this as 

(38) rji(t) < Bi j rji-i(s)ds+pi(t), *£{0,1}. 

Jo 

By monotonicity of 771-i, 

[ rji-i(s)ds < Trji-i(t), 

Jo 

whence 

rjfit) < BiTrji-i(t) + pfit). 

Inserting dm for r]i-i(t) into this, we obtain the following inequalities that hold for every 
t £ [0, T) and i £ {0,1} 

fjfit) < BiBi-iT f rji(s)ds + pi(t) + BiTpi-i(t). 
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Setting 

A = BqBiT, <pi(t) = <pi(t) + BiT!pi-i(t), 

this reads 

(39) rji(t) < A j rji(s)ds + ipi(t), i £ (0,1}. 

Jo 

Defining for i £ (0,1} the absolutely continuous, non-negative, non-decreasing functions 

Hi(t) = f rji(s)ds, 

Jo 

we obtain for almost every t £ [0, T) 

j t Hi(t ) < AHi{t) + cpi(t). 

Since Hi( 0) = 0, we obtain from Lemma [3] for every t € [0, T) and i £ {0,1} that 

Hi(t) < [ e AI ' t ~^ipi(s)ds. 

Jo 

Since, by dTTTTli . (13311 . we have for i £ {0,1}, 

rji(t) < Ai^AHi{t) + tfi(t)') + Bi (AHi-i(t) + $>i_i(i)) + 
we obtain the desired result with the constants 

cl =A t + B Q B 1 Te A '~ iT , 
c]-* =B t + AiBiTe AiT , 
d = e (Ao+Al)T BqBxT, 

f = A, 

9i = A, 

where i £ {0,1}. □ 


3. John’s bounds on the first order quantities and the inverse density of 

CHARACTERISTICS 

3.1. The evolution equations for the first order quantities (cf. [Joh74] '). In this 
subsection, we compute the evolution equations of some first order quantities along an i th 
characteristic. 

So fix i £ {1,..., N}. For k £ {1,..., N}, let (see (I3U1) . (IBTl) J 

d 


k *k 

w = e —u, 
ox 


(40) 

and, for the given i, let 

(41) 

Then, by iRTHl . 

(42) £ w *E( e * fe £ M ) efc = E 

and so, by (121 [I . (HU) . 

jr~ u = Pi ^ w k ek = v l a + pi ^ 


v = e* l -^—u = piW 1 

OZi 


( *fe d 
k 


w k e k , 


d_ 

dzi' 


(43) 

From © and m we have 
d_ 
dsi 


k^i 


d , d d /. . s d 
= A,—u + -u = (\i id -a) —u, 


dx 


( 44 ) 
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whence, by (PI). 
(45) 


dsi 




It then follows from (PI) that 

d 


(46) 

while, using P)) . 


(Jk a ) ek = 'A2 e *^ De ‘ a '>{ e * l ^: u ) ek = - Ai )4 


(47) 


i {iL a ) ek = J2 e * j ( D *' a ')( e * l -^ u ) ek =Y^ w ‘- 


Similarly, from PI) . we get, using 63 , respectively PI) . P)l. 


(48) 

(49) 

(50) 


_d_ 

dsi 

d 


k 

dx ^ = Y. c >^ 


a a .= Pi A wk = ^y+ P i A 


k 

k^i 


(51) 


In the same way, we recover with (S3 from PI) that 

-e = 2_^ \._\ w \ c 3i\ e ^ e j) o e + c iz e )• 


dsi A 7 

k,l 


Finally, we get, using (12U1) . (1T51) . (15U1) . 
(52) 


d _ d ( d \ _ d ( d \ _ <9 _ s—' 

9si V5zi V 9«iV9si V fc P* 2.^1 Cl 


k i i - 
fc W = CjjW + 


^ ^ ('irn/W }pi- 

These are actually not all the first order quantities. As we shall promptly see, ^-io' 
and jf-v * are as well of order one in the derivatives of u. 

OSi 

We assume in the following that K > 2, although equations (1551) . (1571) . (1551) below 
remain valid for K = 1. We start by considering -S-w*. We have from (PI) 


d i 
—w = 

USi 




u + e 


d_( d 
dsi 


But, from CD. and using ©> we have 
d ( d 


-xy 


ds 
whence 


. d 2 d 2 . d 2 d ( d 
ox z axat ox z ox 


JL(y \ 

dsAdx ) 


(£«))■ 

(°S“) = ( A ‘ ld -“)S5“-(s“) 

(yiA 


d \ d 

dx % 


Therefore, by (1511) . as well as (1251) and (1421) . 


_d_ 

dsi 


k,l 

ky6i 


i \ ^ Ai Al i / k / \ i . i k\ \ ^ i k l 

= ) v x w {cu(ek,ei)QW + c kl w ) - 2_^c ki w w 


k,l 


= - C ‘ i (w i ) 2 - (Yfa+y-Y 

A i 

^ A* — Ai 


k 

k^i 


A i — A\ 
A i — Afc 




k,l 

k^i,l^i 


K Ai — Ai 


i \ i k l 

- 1 c fei w w . 
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Thus, defining the coefficients 

(53) 7im = 7im(u) = 7m! 

as in I.Toh74l (notice, however, the difference in the last term in (I54bl) stemming from the 
fact that we used the normalisation condition ||ei||o = 1, whereas John uses |e**| = 1), 

(54a) 

(54b) 

(54c) 


7ii 

i 

7 i 


i _ 1 | i i , \ ^ Xi X m l , \ 

im Tv I ^im > / 7 ^ira j 0 

2 V V Ai - Ai 


'Jim 


1 / A/ Xn 


2 \ Xi — A/ 


. ’4 I 

■Qm H“ 


Am Aj ^ 
Aj Ar 


~Cml 


we can write 
(55) 


5 i i 

= 2^ 7i " 


l m 
W W . 


m i, 

l i,m ^ i, 


Noticing that, by (154cII . 

(56) 7 u = 0, if / + i, 

and using (53), (I54al) . we can write this as 
d 


dsi 


i / i \2 i ^ ^ i l m 
= -Cii(w ) + 2_^ llmW W 


l,m 


(57) 


i(w l f + (X] + XZ 


i l m 
7Im® W . 




We take note that the right-hand side is indeed of first order in derivatives of u, as 
announced. 

It is now easy to compute v We obtain from SD, GD and <E3 


d / i\ 

( 

d \ i 

( 9 

= drM W ) 

= { 

w^r +pA 

Ids™) 

= CiiV l W l + 1 

IE 

i m\ i 

Cim'W I piW 



- Pi{c\i{w X ) 2 ) + 27 * 1 " 


W )W + p, 


( ^2 7 lmW l W m 

l,m 

l*m 


which, by m, immediately simplifies to 


(58) 


d t 
x~ v = 

C/Si 


(E( 27 


tm 1 '-'im 7 


i l m\ 

7 ImW W Pi. 


( E 

m * i l 

l^m 

Again, we take note that the right-hand side is of first order. 

3.2. Some a—priori assumptions on the initial data. Seeking C K -solutions u of the 
Cauchy problem (1 11 1) with initial condition (cf. (1101 ) ) 

(59) u(-,O) = / = 0/ o , 

where fo £ C^’ 1 (R; ( 0)) has support in 7o = [—1,1], we first have to restrict ourselves 

to $ £ (0,1) so that / has range in ( 0). Then, since we assumed that a £ (7 max ^ 2, ' K '} J 
the quantities a, Ai, e* 1 , d, c kl are defined and at least C 1 composed with /. We have 
(cf. |.Toh74h 

/ = O(x)), f = J-f = <D(r)).. 


( 60 ) 
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and thus 

(61) A,(/)-A,(0) = O(tf), e ‘ i (/)-e* , (0) = O(fl), «*(/)-*(0) = 0(0), 
while 

(62) Xi = 0(1), e** — 0(1), ei = 0(l). 

We wish to stress here that we compute lengths of tangent (co-)vectors not with respect 
to the Euclidean norm, but with respect to the appropriate norm induced by the scalar 
product (■, • )o with respect to which the basis {ei(0)}, 6 /i,...,jv} is orthonormal (cf. l[5|)). 
We further obtain that 

*(/)-<&(()) = ©(/) = 0 ( 0 ), 

so that, by (J29D, it suffices to choose 0 small enough so that 

(63) (1 + e)d»(0) < cu(f) < (1 - e)cj»(0) < 0. 

Finally, by (HU1) . (T5U1) and we have for every z£l, 

(64) w\z, 0) = e*V'(z) = O(f') = 0(0). 

Let us introduce now the following quantities: 

(65) Wo = maxsup w‘(z, 0) , Wj~ = maxsup w l (z, 0), 

» z * z 

(66) Wo.o = maxsup |e* I (0)/o(2)| , W 0 + 0 = maxsup(e* I (0)/o(;z)). 

We trivially have 0 < W 0 + < Wo and 0 < Wq 0 < Wo,o- Moreover, (I5H) implies that 

(67) Wo = 0(0). 

Also, we immediately infer from (1621) (using in the first line also (1421) 1 that 

( 68 ) /' = !>*(■> 0)ei(f) = O(W 0 ) 


and 

(69) 


/o = ^(e* l (°)/o)e i (°) = 0(Wo,o). 

i 


Now, since fo is C K and has support in Io = [—1,1], we must have, for all i £ 

(70) 0 = e**(0) fo{z)dzj = j (e*\0)fo(z))dz. 

=o 

Therefore, if fo is non-constant (i.e., if fo is non-zero, since, by continuity, any other 
constant would contradict fo being supported in Io), there must be an index i £ {1,..., N} 
and a point in Io where e**(0)/o is strictly positive. We conclude 

(71) 0 < W+o < W 0 ,o. 


We now establish an explicit lower bound for Wq q. By (1661) . there is an io and a zo such 
that |e* l °(0)/o(zo)| = Wo,o- If e**° (0)/q(zo) > 0, we obtain from (I5S1) that 

(72) W+ 0 >e* i °(0)fl ) (zo) = Wo,o. 


So let us assume that e* l °(0)/o(zo) < 0, and set 


(73) 

where 


to = 


Wo.o 

L 


L = 


sup 


|/o (*0-fo(z" 


(74) 


0 
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is the Lipschitz constant of /q. We then have, following IJoh74l . 

Wo,o = -e* io (0)fo(zo) = -^~ [ e* io (0)ti(z o )dz 

260 J\Z-Z 0 \«L 0 

= y- [ e* l °{0)(fo(z)~ fo(z 0 ))dz+f e* l °(fi)f' 0 (z)dz 

2£ 0 J\z-z 0 K^O 2e ° Jr. „,l>ro 

~ 2 ~ J : e * IO (°)/o(2) rf « 




GSI'llU 27, ,, 

< |rc>l 


i f G3 1 r 

ydy+2 -I Kodz < ^Wo.o + ^w 0 +. 


' Io 


Wo,i 


We conclude that 
(75) 


Wn + n > 


2 L ' 


Since Wo,o < L, for /o(±l) = 0, the right-hand side of (1 75 II is smaller than the right-hand 
side of (1721) . so that (1751) is valid for either sign of e* l ° (O)fb(zo). 

We next show that G3 entails that we can estimate Wo in terms of WqV By USD and 
(PI) . we have for all i 

(76) (e*‘(/) - e^O ))f 0 = (e^fl/o) - e*\0))fb = 0{dW 0 ,o). 

Consequently, choosing d small enough, we can ensure that 

(e**($/o) — e* i (0))/o| < iwbto, Vi G {0,..., N}, 

so that, if we pick i + € (1,..., AC} and z + G Io such that e* I+ (0)fo{z + ) = Wq 0 (cf. (RHUl l. 
we get 

e* i+ (df 0 (z + ))fo(z + ) = W+ 0 +(e* i+ ($f 0 {z + ))-e* i+ { 0))^( Z +) > W+ 0 ~W+ 0 = i W+ 0 . 
But then, since we have WJ 1 " > e** + (/)/' on all Io (cf. (1651) 1. we infer that 


(77) 


W 0 + > e* i+ (df 0 (z + ))(dfb(z+)) > |w 0 + > 0. 


On the other hand, if i and 2 are such that 

(PI), (USD, that 

W 0 = |e* 7 (tf/oP)(tf/o(r))| <0 

= 0(dW 0 ,o + d 2 W 0 ,o) = 0(dW 0 ,o). 
We conclude, using cmi, m , that 


e* l {f{z))f\z) 


= Wo, we recover from PI) . 


e*’(0 )fUz) 


+ d 


(e*“ ( 0 / 0 ( 2 )) — e**(0)^ fo(z) 


-=0( ^ ) =0(1), 

W 0 + \#W+ 0 J 


(78) 


Wo = O(W 0 +). 


3.3. John’s Theorem — the bounds. In this subsection, we recover the bounds ap¬ 
pearing in the proof of the main theorem of F. John ( IJoh741 Section 2]) following closely 
his argumentation. We assume that, for some T > 0, u G C K (R x [0, T]; 73^(0)) is a 
solution of CD, where a and 8 are as in Theorem [5] and •& is small enough so that the 
estimates of Subsection 13.21 hold. 

We start by considering 

(79) W(t) = max sup |tu*(a; 7 , i 7 ) | 

* (x'.t'j 

0<t'<t 
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on [0, T], where w x {x,t) = w l (u(x,t)J is defined by moj. Notice that W is non-decreasing 
in t (it is actually also continuous, since the w/’s are C K ~ 1 and have bounded support in 
x), and that W (0) = Wo. From (15511 . we have for all i £ (1,..., N} and every i£l that 
(80) 


A !«/(*,t)|< 


d i, 
dT™ (M) 


E i l m 

7 ImW W 


l,m 


< r w(tf 


for almost every t £ [0, T], 


where 

r = sup ^ |7i m («)| ■ 

We then apply the following standard Gronwall-type argument (cf., e.g., [MPF9li § XI.13]) 
to obtain that 


(81) 


|iu’(a:,i)| < Y(t) 


for all t 6 [0, T] for which the solution Y(t) of the initial value problem 

. £v = rv 2 

(82) 


f _y = r y * 

< dt 

l V(0) = Wo, 


is defined. Indeed, let Z be a solution of the initial value problem 


d 
dt 

Z( 0) = Wo, 


Z = F(ma x{W,Z}) : 


and assume that there is a time ti for which W(ti) > Z(t\). Then there is an i\ £ 
{1,..., AC} and a z\ £ R, such that 

(83) \w il (X il (z 1 ,t 1 ),t 1 )\=W(t 1 )> Z(t!). 

Let t = infjt | W(s) > Z(s),Vs £ [t, ti]}. Then, necessarily, 

(84) W(t) = Z(t). 

Moreover, on \t, ti], we have ma x{W, Z} = W, so that 

rti 

z{ti) - Z(t) = / r W{s) 2 ds, 


(85) 


z(t)+ / rw(s) 2 ds = z{ti). 


On the other hand, by integrating (I5CT1) for ii along the if 1 characteristic Ci 1 (zi) passing 
through (Xi x (zi,t\),t\ \, we have 

/■T 

w n - w !1 (X il (zi,t),t) < J FW(s) 2 ds. 

This means that, using (1791) . (1841) and (1851) . we get 

/■ti 

\w n (X n {z 1 ,t 1 ),t 1 )\<\w n (X n {z 1 ,t),t}\+ J F W{sfds 

/ t i rt i 

TW(sfds = Z(t) + J VW(s) 2 ds = Z{t!). 

But this contradicts (1551) . and we must have W(t) < Z(t) for all t £ [0,T] for which Z is 
defined. But then Z is actually a solution of (52), which establishes (HD. 

Now, JE2 has the solution 

Wo 


( 86 ) 


Y(t) = 


i - rwot ’ 


t £ 


o, 


u. 

rwoJ 
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Recalling from G9 that the time to after which the characteristic strips 7 Zi are separated 
is 0(1), so that, by (ITTHl . 

rWo = O(W 0 ) = 0(0), 

we see that (El) implies 

(87) \w l (x,t) | < (1 + e)Wo, V® G R, Vt € [0, to], 
provided that 0 is small enough. In other words, using El. 

(88) W(t) = O(W 0 ) = 0(0), Vt G [0, t 0 ]. 

Let us now introduce the following time-dependent bounds defined on [0, T]: 

(89a) l/(t) = max sup | a/r* (ar / , i / ) | , 

0<t'<t 

(89b) S(t) = max sup pi(z'i,s'i), 

i 

0<a' <t 

(89c) J(t) = max sup \v l (z'i, s() , 

i 

Z i eI ° 

0 <s'.<t 

(89d) U(t) = sup |it(:E , ,t , )| . 

o <t'<t 

Like W (t) (see (1791) 1. these quantities are non-decreasing in t (and, actually, also continu¬ 
ous). By (ED, CEO), we see that, initially, 

(90) Pi (zi, 0) = 1, 
and thus, by m, that 

(91) v'(zi,0) ~-w'(zi,0), 

so that we obtain, using that u( ■, 0) = /(•) G B^ 1 ( 0) is supported in Jo, 

(92a) 17(0) = 0, 

(92b) 5(0) = 1, 

(92c) J(0) = Wo, 

(92d) U( 0) < 5, and, by I®, U(0) = 0(0). 

Next, we want to estimate these quantities on [0, Jo]■ For V, we make use of (1801) and 
the fact that, for any (x',t') qL IZi the corresponding characteristic coordinates (z(,s() = 
(«'(*', t'), s' (x', t')) fulfil z'i 0 Jo so that w‘(z', 0) = 0. Then we have, using lUHl) . (1551) . 

(93) V(t) = 0(WS), Vt G [0, t 0 ]- 

We then consider pi. By the second to last equality in EJ- we see that 

(94) pi(zi,Si) = 0(W(si)pi(zi,Si )), 

from which we immediately conclude that, for as long as u is defined and C K (i.e., for 
t G [0, T] according to our assumption), so that W is bounded, we have, using the standard 
Gronwall-inequality (Lemma [3J , 

pi(zi,t) > pi(zi,0) exp^-0 (tw (t))j , 

=i 

and hence 

(95) Pi(-,t)> 0, V* G (1,..., JV}, Vt G [0,T\. 

In particular, we have for t G [0, to], using (TT51) . (1551) . 

pi(zi,t) > exp [-0(0)), 
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so that we get for small enough d, 

(96) pi{ • ,t) > 1 — e > 0, Vt £ [0, to]. 

Similarly, (1551) implies that pi(zi,t) < exp ^C>(tW(t)^ , so that, by (1551) . (I89bl) . 

(97) S(t) < (1 + e), Vt £ [0, t 0 ], 
for sufficiently small d. In other words, 

(98) 5(1) =0(1), Vt E [0, t 0 ]. 

Regarding v', we obtain from (1551) . using ED , (1751) . (USED , that 

-^v\zi,Si) = 0(piW(t) 2 ) = 0(S(t)W(t) 2 ), 

which, by ifisp. ED, ED, EiiD, ED, implies after integrating and taking supremums 
that 

(99) J(t) = J( 0) + 0(Wo) = W 0 + O(dW 0 ) = O(W 0 ), Vt € [0, t 0 ]. 

Finally, since 

px 

u{x,t) = / 


!x n (- l,t) 


-^-u(x',t)dx', 


we obtain from ED, 


( 100 ) 


i(x,t)\ = ol ^2 


r x 1 (i,t) 


|w l (a; , ,t)| dx 


>x N (- i,t) 


Then, by (I89dl) . (1T51) . (1751) . (1T51) . 

(101) l7(t) = 0^(2+(Ai(0)-Ajv(0))t)w(t)j = O{W 0 )=O(d), Vt£[0,t o ], 


so that 


(102) U(t) < S, Vt E [0, t 0 ], 
for sufficiently small d. 

Using (1571) . we conclude from (1551) . (1551) and (1571) . and (1551) . respectively, that V, S, 
and J stay close to their initial values ED , as long as t E [0, to] and d is suitably small. 
Also, under those assumptions, the solution u remains in the Al-ball B™( 0) where it lay 
initially. 

We next consider times t E [to,T]. Let ( x,t ) E IZi for some i be given. Then the last 
equality in (1521) gives along the i th characteristic Ci(zi) through ( x,t ) that 

(103) -^-pi(zi,Si) = £)(j(si)+ V(si)S(si)), Vsi E [to,T], 

OSi 

where we have used that, by JT5J, (jZm f) IZi) f'|{(a:', t') | t' E [to,T]} = 0, Vm ^ i. 
Integrating, and using ED, we obtain 

Pi( Zi ,t) = 0(1 +tJ{t) +tV(t)S{t)), 

so that, by varying i and x such that ( x,t ) G IZi ('[{(a/,!') | t' E [to,? 1 ]}, we have 

(104) S{t) = 0(l + tJ{t)+tV{t)S(t)), Vt E [0,T], 

which is compatible with ED - Next, we argue in the same way for v z . Fixing again 
(x,t) E IZi for some i and using (1551) . we get 

(105) -^-v\zi,Si) — 0{V(si)J(si) + V(si) 2 S(si)) , Vs, E [to, t]- 

Integrating, using EJ, and varying i and x, we arrive at 

(106) J(t) =O(W 0 + tV(t)J{t)+tV(t) 2 S(t)), 


Vt E [0, T\. 
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From (tnm and (ITM1) . we see that it will be fundamental to estimate V(t). So fix i £ 
{1,...,7V} and assume that ( x,t ) ^ IZi with 7 £ [to,T], (1551) then gives along the i th 
characteristic Ci(zi) passing through (x,t), 

(107) ^-W ! (Xi(2i,Si),Si) = ^(7 ImW'W™) Si), Si) . 

I ,m 

For the given Zi, let 

( 108 ) n£*(7) = jV £ [ 0 , 7 ] en k y ke{i,...,N}. 

Clearly, FI?*(7) = 0, since Zi ^ Io. Also, by definition (fTHl) of to, so that (fTTJl) holds, 
(n^(7) P| IIf4(7)) P[to, T 1 ] = 0. Therefore, using Zi Io, whence w l (zi, 0) = 0, as well as 
(1551) . (I89al) . (1551) . we obtain after integration along Ci(zi) 

(109) \w x (zi,t)\ = o(wq +tV(t) 2 + j \w k (Xi(zi,t'),t') | dt' j . 

' nk ' (t) ' 

Observe that 

® In =0, if Mt, 

was crucial for this estimate. Now, for each t' £ 11^(7), we have (Xi(zi,t'),t) £ 7 Z k . 
Parametrising Ci(zi) {~\!Z k by y k £ 7^(7) for some 7^(7) C Io (recall Subsection 12.41) . we 
get from (I22bl) . using that yt = zi is constant along Ci(zi), 

Pk(yk,t l (y i ,y k )) 


(HO) 


\w k (Xi(zi,t'),t') | dt' = 


|Ai Afe| 


w k (yk,t'(yi,y k )) 


dy k 


\v k {yk,t'{ yi ,y k ))\ 

-p;— ^- dy t =0(J( .». 


Consequently, we get from GSB) , (I55al) after varying i and x such that ( x,t ) ^ IZi, 

(111) V (7) = 0(Wq + t.V(t) 2 + V(7) J(7)), 7 £ [0, T\. 

Finally, going back to (1 1001) and sub-dividing the x-interval [Xat(—1,7), Xi (1,7)] for 
7 £ [7o,T] into regions where ( x,t ) £ 7 Zi for some i and regions where (x,t) 0 (J fc 1Z k , we 
obtain with mu 

\u(x, 7) | = 0^(X 1 {l,t)-X N {-l,t))V{t)+J(t)^j =C^2+(Ai(0)-Ajv(0))7)v(7)+J(7)^, 

where we have used that, by m, changing variables from ( x',t ) to ( 2 ',7) yields for 
constant 7 

r x i( i,‘) r 1 r 

\w\x ,t)\dx = / \w l (Xi(z'i, 7), 7) | pi(z'i, t)dz'i = / \v\z'i,t)\dz'i = 0[j(t)). 




'Io 


Consequently, we obtain taking supremums 

(112) U(t) = O(V (7) + tV (7) + 7(7)), 7 £ [0, T], 

where we have used that, for 7 £ [0,7o], (11011) . (ITiJl) . (TUI) . (IM1) and (189cl) allow us to 
estimate U{t) by J(7). 

Summarising, we have from cm. tm cm) and cm , omitting the explicit 7-de¬ 
pendence, 

(113a) V = 0(Wo +tV -V + VJ), 

(113b) S = 0(1 +tJ + tV ■ S), 

(113c) J = O(W 0 + tV ■ J + tV-VS), 

(113d) U = 0(V + tV + J), 


which are valid for all C^-solutions u of dill) for which u £ 77^(0) for all 7 £ [0, T], 
provided •& is small enough. 
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We now perform a bootstrap argument. We let t increase from zero so as to keep the 
following inequalities satisfied: 


(114a) 

(114b) 


tv < J < yfd, V <■&, 


u <V&. 


When t = 0, the estimates (111 4b ft hold by (19211 . (1671) . for all sufficiently small $. Notice also 
that the right-hand side of (I114al) is a constant (depending only on a and e) which is finite, 
since we assumed the system to be genuinely non linear. We then get, for sufficiently small 
0 , 

V = 0(Wg + \f&V + VtiV) => V = 0(Wq), 

S = 0(1 + tJ + VtfS) ==> S = 0(l+tJ), 

J = O(w 0 + VtfJ + VtiVS) =4> J = O(W 0 + VtfVS), 

so that 

tv = o(twS), 

VS = 0(V + tv ■ J) = 0(V + yf&J). 

Then, using (I114all . (1751) . (1571) . we obtain 

V = O(W 0 2 ) = 0(tf 2 ), 

tv = 0(tW$ ■ Wo) = O(W 0 ) = 0(0), 

J = O(w 0 + Viiv + tiJ) => J = O(W 0 + VtiWo) = O(W 0 ) = 0(i9). 

Finally, this gives 

U = 0(Wo + Wo + Wo) = O(W 0 ) = 0{d), 

so that, indeed, u stays in B^( 0), provided t) is small enough. So we have established 
that, as long as t is so small as to keep the conditions cnn fulfilled, we have the stronger 
estimates 


(115) tV = O(W 0 ) J = O(W 0 ) V = 0(W 2 ) U = O(W 0 ) 

= OW, =o(t>), =0( tf 2 ), =<3(tf), 

which, therefore, hold as long as (I114al) does, the assumptions (I114bl> being unnecessary. 
For later reference, we also list the following additional estimates that thus hold (using, 
as well, G2D and my- 

S = 0( 1 + tJ) = 0(1 + tWo) = 0(1 + tw 0 +) = 0(1), 

( 116 ) o o 

VS = 0(V + tV ■ J) = 0(W$ + ■dWo) = 0(dWo) = 0(tf 2 ). 


and (HTB1) . together with (1351) . establish statement 0 of Theorem [2] (with the excep¬ 
tion of the bounds on pi and v l outside IZi, which are obtained in the next subsection). 

We now focus our attention on the inverse density pi of the i th characteristics, thereby 
also establishing the second part of F. John’s theorem n.loh741 Section 2]). 


3.4. The inverse density of characteristics. In this subsection, we analyse more in 
detail the behaviour of the inverse density of characteristics pi. We assume that u £ 
C k (R x [0,t*);Bf(0)) is a solution of the strictly hyperbolic and genuinely non-linear 
system © with initial condition 00 satisfying the additional assumptions set out in 
Subsection where t* is the maximal time of existence of u. Among other things, we 
will give an explicit estimate for the size of t*, as well as establish that pi has to vanish 
somewhere at that time. This then recovers the second part of F. John’s main theorem 
f 1 Joh74l Section 2]), since ( 14111 . the bound on J in ( 11151 ) and the lower bound on v l in 
cc© below imply that w l has to blow up as t approaches t*. 
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Our starting point is the evolution equation <E3i of the inverse density pi along an * th 
characteristic, 


(I521l 


d i i 

J—Pi = CuV + 
OSi 


X 


i m \ 
Cirri'^ J Pit 


which, given our sign convention (I2il (i cli < 0), we rewrite as 

d 


(117) 


ds 


■Pi = -\cli\v + 


X 


i m \ 

('im'M ) Pi’ 


We want to show, in order to obtain Theorem [21 (lull , that there is an index i + £ {1,..., AC} 
and a point z + £ Iq = [— 1,1] such that p^+ (z + ,tj tends to zero as t /* t*, and that 

( 118 ) ue[n,%\, 

where 

1 ^ 1 


(119) To = min ■ 


and To = max ■ 


(1 + £) 3 I4(0)|W+ c (i- e ) 4 |4(0)|w+- 

Notice that tW$ then fulfils (llldall for all t £ [0,To], so that the estimates (1X151) . (ITTB 
hold for any t £ [0, T], T < To, for which we have a C K solution tionRx [0, T]. 

First, let us define 


(120) a(t) = max sup 

«'e/ 0 m * 

0 <s'<t 


y, C im (u(Xi(z', ( Xi(z ', s'), s') 


so that, by (11171) . 


-\cu\v - api < — pi < -\cu\v + api. 


Then, if v l > 0, we have by (O for small enough d, 

(121) - (1 + e)|4(0)|r/ - api < T^pi < -(1 - e)|4(0)X + api. 

Now, by (l88l) . (I89al) . (11151) . (1671) . we have for all t £ [0, min{t*,To}), 

[ a{s)ds = O(W 0 + tV) = 0{d), 


so that, for sufficiently small d, 

exp / a(s)ds I < (1 + e), implying exp — / a(s)ds I > (1 — e). 


Consequently, by the standard Gronwall-inequality (Lemma [3j), (1901) and (11211) imply for 
sufficiently small d, and for as long as v l > 0 along the i th characteristic starting at Zi, 
that 

(122) (1 - e) (l - (1 + e) 2 |c«(0)| [ v\zi,s)ds] < pi(zi,t) 


< (1 + e) l-(l-e) 2 |4(0)| / v\ Zi ,s)ds\. 


In particular, 

(123) Pi < (1 + e) 

as long as v l > 0, provided d is small enough. 

Now recall (1551) . 

d 


((551) 


dsi 


(E( 2 ^ 


4774 I '-'4771 / 


m^i 


( E 

l, m 
l^m 


w l w m I pi. 
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Letting 
(124) 

Pi (t ) = max sup 
*'e/ 0 

0 <s'<t 


P 2 (t)= max sup 

0<s / <t l^m 

we have, whenever u 1 > 0, and using dm , 

& 


Y. (27 im +dm) («(-Xi(z^sO,s'))t(; m (.Yi(z\s'),s') 
y 7?m s'), s')) w' (Xi(«, s'), s') w m (Xi(z , s'), s') 


(125) 


9s, 


u* > -^lU 1 - /?2Pi > ~Plv' - (1 + e)/?2. 


So, again by Gronwall and arguing as before, taking into account dH3 , dm , dS3> , we get 
for sufficiently small d that 

v'(zi,t) > (1 - e) (w\zi, 0) - (1 + e) 2 [ p 2 (s)ds], 


provided that v' l (zi,s) > 0, Vs G [0,t]. In particular, taking i + G {1,..., IV} and z + G Io 
such that v I+ ( 2 + , 0 ) = w* + (z + ,0) = W 0 + > 0 (recall d2Bll . dTTTll . dSSJl), we obtain 

v l+ (z + ,t) > (1 - e) ^Wq + - (1 + e) 2 J p 2 (s)ds S j . 

Observing that, by dED, dm, dszD, dza, 

[ p 2 (s)ds = 0 (Wq + tV 2 ) = 0(Wo + dWo) = 0(dW£), 

Jo 

we obtain, for sufficiently small 

(126) v i+ (z+,t) > (1 -e)(W+-eW+) > (l-e) 2 W+. 

In particular, v' l+ (z + ,t) > 0 for all t G [0, min{t*, To}), and we did not need to assume 
v 1+ > 0 along Ci+(z + ). Also, going back to dm , we have 

(127) p i+ {z + ,t) < (l + e)^l-(l-e) 4 |c^ i+ (0)|tW / o + ^, Vt G [0, min{f*, T 0 }), 


from which it is immediate that lim t/'T + Pi+{z + ,t) = 0 for some 


(128) 


T i+ < 


(1 - eY 




-(Q)Fo + 


< Tq. 


Since we already argued at the beginning of this subsection that u cannot exist after the 
time T t + when p t + becomes zero, this establishes the upper bound for t* in (fTTSl) . which 
is similar to the one obtained by John. 

To obtain the lower bound on f*, we first consider again dSSD, from which we immedi¬ 
ately obtain by (141 ll . (tHHl) . (IM1) . (11811 . (IHH . (11151) . (IllfiD . (RT7T) 

v i (z i ,t)-w i (z i , 0) < 0(Wq +tV J + tV 2 S) = O(dW 0 + dW 0 + d 2 W 0 ) = O(dW 0 ), 


i.e., 

(129) v i (z i ,t)<w i {z i ,0) + 0(dWo). 

Similarly, and for later reference, we obtain 

(130) v\ Zi ,t) > w\zi, 0) -O(dW 0 ). 
Consequently, we can assume by dssD, m~ that, for suitably small d, 

(131) v i <(l+e)W 0 +, Vi G {1,..., Nj. 

Now, going back to dm, we have with dm 

d { i 

"o Pi > -\ca\(v )+ - api, 

OSi 
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for v l < (v t )+, whence —v l > — (v t )+, where (?/)+ denotes the positive part of v'. Thus, 
by (1631) and for •& small enough, 


(132) 


-S-Pi > -(1 + e)|c«(0)|(u l )+ - api. 

USi 


Then, by Gronwall (Lemma [3]), using (1881) . (189aD . (I115D . (1901) . we get for sufficiently small 
•& and all t G [0,f*) (we already showed that t* < To) 

(133) pi(zi,t) > (l-e)^l-(l + e) 2 |c)i(0)|^ {y\zi, s)) + ds^ . 

Thus, (11311) implies that 

Pi(-,t) > (l-e)(l-(l + e) 3 \cU0)\tW+). 

Therefore, 

pi(-,t) > 0 as long as t< 


(1 + e) 3 |cb(0)|W / g"' 

We conclude that the inverse density of characteristics in each characteristic strip remains 
bounded away from zero as long as t < To (cf. (11191) 1. Thus, from (1421) . (1411) . (11151) . we 
see that u G (0) and remains bounded, for all times in [0, min{f*, To}) • Thus, if 
t* < To, we should be able to extend u across t = t*, which would contradict the definition 
of t* as the maximal time of existence of the solution u. Therefore, also the lower bound 
in (fITHI) must hold. This finishes the proof of Theorem El (HI)) . 

We next wish to show that pi remains bounded away from zero along every j th char¬ 
acteristic ending (that is, for t approaching t*) in a point outside 7 Zi. Before doing so, 
however, we quickly analyse the behaviour of pi outside 7 Zi. From (1521) . we have for every 
Zi To 

-a{zi,t)pi(zi,t) < —pi{zi,t) < a(zi,t)pi(zi,t), 


dsi' 


where 


a(zi,t) = sup 

o <s'<t 


Cim (u(Xi(Zi,s'),s')^W m ( Xi(Zi, s'), s') 


which is non-decreasing in t. Then, by (f7<71) . (IH51) . (IHT71) . (11081) . ( 11101 ) . (1(71) . (11151) . we get 
for all Zi qL Io, t G [0, t*), 

[ a(zi,s)ds = ofwo+tV + ^2 [ \w k (Xi(zi,t'),t') | dt' 

Jo V fc 

= o(w 0 + tv + j) = o(P + 1 ? + d) = o( d). 


Thus, using Gronwall with ([HD , we can ensure that, if & is chosen suitably small, we have 


(134) 1 — e < Pi(zi,t) < 1 + e, Mzi qL /o, Vt G [0, t*). 
We observe that this then also implies that, using OH), (UT5l) . 

(135) = O(W 0 2 ), Vzi #I 0 ,Vt€ [0,u). 


This thus establishes the remaining part of Theorem [5] Q). 

Let now (x, t) G (lx [0,t*)j \(J fc 7 Zk and consider the j th characteristic Cj(zj), Zj ^ Io, 
passing through (x, t). Assume further that i G {1,..., TV} is such that Cj(zj) P| IZi 7 ^ 0. 
Fixing re G ^0, ^ (for the sake of definiteness, let re = ^z£l), we face two possibilities. 

Either pi > re all along Cj(zj) (~}!Zi, or there is a point on Cj(z-j) f)!Zi where pi < re. 
Having nothing to prove in the first case, we focus our attention on the second case and 
show that pi cannot become too small along Cj(zj) f)IZi. Recalling Subsection 12.41 let us 
parametrise Cj(zj) {~\!Zi by ip G Io (( 27 , 0 ) and (x,t) both lie outside IZi), so that 


Cj{zj) IZi — (yi,yj)),t ( yi,yj )^ 


Vi = Zj,Vi € Io 


By our assumption, there must be a yt £ Io such that pi(yi,t'(yt, Zj)) < re, so let us study 
pi along the * th characteristic Ci(yi) starting at yt. Our first claim is that pi must be 
strictly decreasing whenever it is smaller than re, at least when 1 ? is sufficiently small. To 
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see this, we first establish that v 1 needs to enjoy a positive initial lower bound in this case. 
For assume that, for any Zi £ To, 

(136) 'S’ ^-dl±|l(i-2 e -«)vr„ + 

(1+e) 2 max;|c‘ ; (0)| T 0 (1 - e) 

Then, by (11291) . (1751) . 

v\ Zi ,Si) <w\ Zi ,Q)+O(VW 0 ) < A ■ 7^4(1 -e~K)W+, Vsi £ [0,t*), 

To (1- — e ) 

if i? is chosen to be suitably small (but independently of k). Since this bound holds also 
for the positive part (w l )+ of v\ GSB then implies that 

Pi{zi,Si) > (l-e)fl- (1 + e) 2 |c' i (0)|-=- • + £ | (1 - e - K)giW 0 + Y 


To (1 - e) 

which, since Si < f* < To, gives by (fTlTfl) 

(1 + £ ) 3 | c ii(0)|(l — e — k)Wq 


Pi{zi,Si) > (1-e) 1- 


>(l-e) 

We conclude that 


(1 + e) 3 (max(|4(0)|) W 0 + (l - e) 

1 — e — (1 — e — k) 

1 — e 


Vsi £ [0,t*). 


w\yi, 0) > -=- • (i _ 2e - /s)W+ 


To (1 - e) 

so that pi < k is possible along C;(j/i). But then, (11301) and (1781) imply that, for sufficiently 
(but independently of k) small •&, 

v'(yi,Si ) > -=-77^-4(1 - 3e - k)W 0 + , Vsi £ [0,i*). 

Jo l 1 ~ E ) 

We thus see that, by (fimi and using (1551) . (I89al) . we have 

~-S~Pi(yi,Si) > (l + e)|cF(0)|^(l-3e- K )W 0 + - K O(l/), 

OSi J Q 

whenever pi < k. By (11151) . V = 0(Wq), hence, using (1571) . (1751) . we can estimate the 
second term on the right-hand side for suitably small •& (again, independently of k) by a 
small fraction of the first term, say by 

-(l + E)|d 1 (0)|-=-e K W 0 + , 
j- 0 

so that we obtain 

(137) - pi(yi,Si ) > (1 +e)|cii(0)|-p-(l - 3e - (1 + s)k)W^, 

OSi j[ Q 

which is valid for all Si £ [0, t*) for which pi{yi,si) < k. This means that, once on an i th 
characteristic pi drops below k, it thereafter steadily decreases. It follows that 

pi(yi,Si)<K and (11371) holds. for all s; £ [t'(yi, yj), t*). 

As a consequence, letting AT = t — t'(yi, Zj) > 0, we have 

rp 

Pi{yi,t- AT) > pi(yi,t) +(1 + e)|cii(0)|-=-(l - 3e - (1 + e)n) W / 0 + AT 


(138) 


rp 

> (1 + e)|cj i (0)|-=-(l - 3e - (1 + e)k)W^AT. 
To 


Therefore, if we can find a lower bound on AT, we will have established that pi remains 
bounded away from zero along Cj(zj) {^\7ti- To do so, we will make use of the fact that 
the seed /o for the initial data / is C 1 ’ 1 and compactly supported in Iq = [—1,1], so that 
fo and /o vanish at the boundary dlo = {±1} of Iq. For simplicity, we restrict ourselves 
to the case i > j, the case i < j being shown in an entirely similar manner, requiring 
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only a few obvious modifications. Under this additional assumption, we must have for our 
initially picked point ( x,t ) that x > Xi(l,t), i.e., ( x,t ) lies on the right of IZi. We then 
denote by 

( 139 ) di(zi,Si) = Xi(l, Si) - Xi(zi,Si), Zi 6 I 0 ,Si € [ 0 , U), 

the ^-distance at time Si of the i th characteristic Ci(zi) emanating from Zi £ To to the 
right boundary Xi(l,t) of the i th strip IZi. Noticing that (161H implies for sufficiently small 
1 9 that 

!-e = (! - e )|| e *‘(°)|| 0 - || e* l (/)||o ^ (! + e )|| e * l (°)|| 0 = 1 + e > 

we see from m . dsgii . the fact that fo is supported in 7o, and (El, that 


p («i,0) = 


< (l + e)||/'(^)|| 0 


< (1 + e)i?||/o(«i)|| 0 = (l+e)tf||/o(2i) - /o(l)|| 0 


< (1 +e)'9L\zi — 1| = (1 + (zi, 0). 

As a consequence, since (HMD implies that pi > k, we have: 


if d r i(zi, 0) < 


To (1 - 2e - k) W+ 


To (1 - e) #L 
We note that, in particular, we must have 


then pi(zi, Si) > k, Vs, e [0,i*). 


To (1-2 e-k)W+ 

W)> tT d- e) it- 


Moreover, applying the mean value theorem to the definition (USD of pi, we infer that: 

then di(zi,Si) > Kdl(zi, 0), Vsj£[0,t*). 


if d r (z .Q)^ ^ (1 - 2g - ^) 

d t [ Zl ,U)< _ (1 _ e) 


Since characteristics of the same family never cross (cf. (11411 1. we conclude that 
(140) 

By lfT51) . ©, we have 


jrr w To (1-2 £-k)W+ w rn . 

di (yi, Si) > -=-(1^7)— ~&l7 K ’ Vsi G I0,7*). 


x - Xi(yi,t'{yi,Zj)) < XjAT 


and 


Xi{yi,t) - Xi(yi,t'{yi,Zj )) > A^A T, 
so that (recall i > j, so that, by ©, A, < A, < A. < A j) 


AT ■> x-Xj{yi,t) > d r i(yi,t) 


A,-A, 


A,-A, 


Setting then 
(141) 

we recover from (HU that 


Yl — Ai AiV’ 


AT >^. ^Z 2 £ - K) ^Ln. 


To 


Inserting this into HMD, we get 
Pifat- AT)>|t±|||4(0)| 


To 

Tb 




(1 — e) 0E L 

((1 — 3e) — (1 +s)k)( 1 — 2e — k)k 


i>TL 




i?L ’ 
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where the last inequality follows from our choice of k = (recall also that e was 

supposed to be very small here, e < 1/10 would have been enough). Since yi € Io was 
arbitrary (but such that pi(xji,t' (yi, Zj)) < k), we finally obtain using dm), El. <nm 


(142) pi > min | ^ ^ , dC', iE j> along Cj (z IZi, if i? is small enough, 


where 


1 (1 g ) 10 |4(0)l ( m in fc |ci; fc (0)| \ 2 (Wo,o) 4 

432 (1 + e) 5 S ^ maxj |cj,(0)| ) L 3 


is a constant depending only on i. e, a, 5 and the C7 1 ’ 1 —norm of /q. 


4. The bounds on the second order quantities 

4.1. The evolution equations for the second order quantities. We use (1571) . (15211 
and (15811 to obtain the evolution equations along i th characteristics for -§^w x , pi = p, 
and v l = respectively, obtaining second order quantities on the right-hand side. 

However, for pi and v l , these will depend on = -^ppj and ^u 3 = -g^-v 3 ( j £ 
{1, ..., W} \ {i}). Recall from Subsection 12.41 that yi is held fixed when taking in the 

dehnitions E} of ^Tj and ( V, so that the latter are the derivatives along i th charac¬ 
teristics of pj and v 3 , respectively, in bi-characteristic coordinates. As we will need to 
bound and ^uj 3 inside IZi, we compute their evolution equations in bi-characteristic 
coordinates along j th characteristics. Using those, we will be able to obtain bounds inside 
7 Zi that depend on bounds at the boundary of IZi. 

First, remark that, using EJ, (HU, we have 

^ _9__9_ d d __ ( d \ d _ /y^ j k \ d_ 

dsidx dx dsi Vcte V dx y2_^, c '‘ kU ) Q x ' 

k 


(143) 


. dsi ’ 


JL 

dx. 


Therefore, El gives, taking into account (El and (1531) (recall also the notation El from 

12.5B , 


i i ( & i\ (\ ^ i m\ ( & i\ 

= -c liW 


dx 


Ci(w') 2 +(J2 27L 


w )w" + 


Y ^ 


l,m 

l^m 
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i i ( & i\ f \ ^ i m\ f & i\ 

- CiiW fe W ME" )\d-x W ) 

m 

m^i 

- ( D eiCii){w 1 ) 3 - (YjD^n?) ( w I ) 2 - 2cuw l (J^w 1 ') 


p^i 


+ ( Y { 2D ^irn)w Tn \w' 1 ) 2 + ( Y, ( 2 - D «= p 7im)w m W P ^W l 


n.^i,p^i 


+ (E 2 ^{iL wm )) wi + (E 2 ^ m ) (lr l ) 


+ +^y^X De rnllrn)w\w rn f 


l^m 


l , m 
l^m 


+Y( De p ’ ylim ) wlw " lwP + 


l,m,p 

l^m,p^l,p^rr, 


2 E^(^ wm )’ 


l^m 


(144) 

d f d A 
d Si \dx W ) 


- 3c >’(^ W ‘) - (y cL« m ) (^) + (E 2 ^ m ) 

rrlji rn=fii 

+(^ 2 M^ wm )) wi+ 2 ^ Lw iL wrn ) 

l^m 

- (D ei Cu) (w') 3 — ^^(/?e p C«)w P ^ (W 1 ) 2 + (£ (2De i 7im)w mS )(w I ) 




+ (^(2-De p 7im)w m W P ^W I + (Y( D ^i'yim) wl 'W m 

m >P Z ,m 

l^m 

+ ^(Z3e p 7!m)w‘w m W P + 2 ^ (-De m 7? m ) W (w™ ) 2 . 


l,m,p 

l^i,m^i,p^i 

l^m,p^l,py6rr, 


l^m 


We will make use of this equation to show the boundedness of -§^w x (iG {1,..., N}) in 
the region (R x [0, £*)J i- e -> outside a ll the characteristic strips, using (1142[) . As 

will become apparent in the next subsection, we should also control -§^pi (i £ {1,..., N}) 
outside all the characteristic strips. But this will follow from the bounds on -§-w x and the 
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following evolution equation of -§^pi along an i th characteristic. 

<9 / <9 

dsi V dx 


\ CM} 

( 9 

A,)( 

( d \ 


< d \ 


) = - 

\dx 

K&c pi ) 


ydJ~i Pi ) 


mm 

( 9 

\ / 

' d , A 

, d , 

^ ,\EU 

d 


\dx 

Pi )\ 

,dx Xl J 

1 + ad 

Kdzi Xi ) 

dz, 

nu 

= pi 

~( 

dx V 

X 

m 

i m 

C-im'W 

n| 

i ^ ^ {^D e i Cim 
l,m 

)w’ 


(145) 


E i ( d m \ 
Cim \dx W ) 

m l,m m 

rr^i 

+ [D^Cu )^ 1 ) 2 Pi+(^ { D ^m C ii + De i C irn )w rn ^W l pi 

+ ^(De m cl m )(^ m ) 2 j/)i+ ( ^2 (D ei C irn ) wlwrn ^Pi‘ 

rr^i 

l^m 

Before computing the evolution equations for yi and v', we want to point out the 
following useful consequence of eqs. (I22al) : 

<9 <9 pj <9 

dzj dyj A i — Xj dsj ’ 


(146) 


j6{l,...,JV}\{i}. 


We can use this to express pj and v 3 in terms of quantities we will be able to control. We 
have along an * th characteristic and for all j ^ i, 


(147a) 

(147b) 


d 


d 


Pj - ur pj ~ ~ pj 1 


dyj 


V = — - V J = — - V J + 


fhrii j«)" <v + a^x(sj«)’ 

£L _ ( j. Pi (A 


dz 


dyj 


= y + 


■(&4 


Also, note that we can write for every j £ {1,..., N} \ {*}, 

Pi f d 


d 


d 


(148) 


j w 3 Pi ( ^ 

= pi — W J = — ' 

OZi OX Pj 

= 4 (-w j pj + u j ). 

Pj 




Using then (fl47l) . and taking into account m, 63, issi)- we obtain 


_d_ 

dzi 


(149) 


W J 1 

% + 

1 ( 


II 

"O | 

'o.lo| 

— Pi 

Pj 

‘V 

W 3 ( 

‘4 + 

1 ( 

.. . 

A pl 

—Pi 

Pj 

l) bJ 3 

W 3 ( 

% + 

1 ( 

.. . 


— Pi 

Pj 

'hr 7 

W 7 ( 

l) Tj + 

1 ( 

, 

A pl 

— Pi 
Pj 


W 3 ( 

l) Tj + 

1 ( 


“ ”p? Pi 

—Pi 

Pj 



Pi 


Pj(^j ^i) 




d 


Xj X 


A j Xi 




'y 3 pq w p w q 


Hi (w 1 ) 2 + 

=0 


(E 2 ^ 


p#i 


3 W P )w l + ^ 7 l q W P 

P,Q 


^ 3 


rx (E 2 < wP ) vi + n (E ^ 


P,Q 

p^i,qy6i 


Now, regarding the evolution equation of pi , we get from dS2D, taking into account 


(TS1) . 


A, 


= JL, 

( d \ 

dsi 

K. dzi Pl ) 

dzi 

Kdsi Pl ) 
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x m , 

, Cirri'^ ) Pi 


m^i 

= {DeiCii )^') 2 + (y^(D ep C ii )w P ') pjV +Cu(-^ r V lS ] 
v 1 

V ^i 

+ {j2^Am)w m y i V i +(^ J2 {De p cL)w™W P y 2 

m m,p 

m^i m^i,p^i 

+ (E ^{-k wrn ))p i + (E (^) 

m^i m^i 

= (E + c\jV l + ( D ei c]i)(v *) 2 + ( y (-De p c‘ m )«; m «; p jp? 

m m,p 

m^i,p^i 

+ I ((-De m 4) + (•D ei 4»)) u,m ^P iV ’ + (E c -(^- w i)^- 

m^i 


Using (fUM to re-express the last term, we obtain 


(150) 

_ 9 _ 

dsi 




Pi = Cim® JPi + <ki v - 

"4 


(E 


-ML 


Tm I pi + 


/ £i2L 

Pm 


(i) m\ 2 

^ Pi 


+ (D e< 4)K) 2 + (^(^cD^^P? + (E X?Ta-(E^> V ) W 

m,p \ m 171 1 p,q 

m^i,p^i m^i p^i,q^i 

]T ( (D em 4) + (D ei cL)) Pi^ + X~~Xi (E^ 


m p \ . 

w y I I p*i> . 


p^i 


Similarly, we get from l|58l). taking into account (S3]), (ESI), 


9 i _ _d_(_d_ i\ _ _d_/_d_ A 
dsi V dsi \dzi V ) dzi V ds4 / 


im* 


■ * \ m \ i , ( \ ^ i l m 

zm+Cim)W )U + I ^ 7 Im^ W 

l,m 
l^m 


)pi 



\ rn \ 

( d i\ 

)w j 

\dzi ) 


l^m 


l^i,m^i,p^i 

l^m 


+ 27L4 ) Pi + (E74^4-) ( Ap.) 


l,m 

l^m 


l,n 
l^m 
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(£ 


l^m 


i l m 
7 lmU> W 


Pi + ( 27 im + 4n)w m )l/ 


n^i 


+ (?D eil \m + ^ ei C- m )w; m ^ (v*) 2 + f^ (De p ^lm)'W l W rn W P ^) pj 


n^i 


(2> 


l,m,p 
l^i,m^i,p^i 
l*m 


eiTmp + 2£>e p 7im + 


i \ m p\ 
im)W W P \ 


piV 




+ ^E 27 imW 1 ^ P* + ^E ( 27 * m + C * m ) ( 


' z, 

l^m 

Using to re-express the last two terms, we obtain 
(151) 
d_ 
dsi 


dzi 


w w . 


fyU lmW l W rn )fM + (^ ( 27 ** m + C i^) U,m ) !/I 
l,m m 


l^m 


l,m 

l^m 


E 


(2 7 | m + cL)w m (i) 


3 )p^+ (E 


pl 

Cii 




Id pjt; 


+ (2/? ei7 L + A iC L> m ) (^) 2 + (E 2 ]33 a 7 (23 27 ‘> p )) («*) 


+ (E^-H"^ + (E £^-(E 7 ^ 


l,m,p 
l^i,m^i,p^i 
l^m 


% „,,Z 


p^i 


qU) U) Pi 


l,n 
l^m 


p,q 

p^i,q^i 


+ {YXDej^ + 2A p 7im + D ep c‘ im )w rn w p J PiV 1 

m,p 

m^i,p^i 

ST' 27 * m + c( m / 

“i Am — Ai -i 


+ 




m p q \ . 

7 pq W^W H I I PiV . 


l*m 


n^i 


p,q 


We have thus obtained a inhomogeneous system for (p,i,Vi) whose coefficients depend 
only on first order quantities, except for some terms in the inhomogeneous part, which 
also depend on the second order quantities ( 7 Vj, 7 \iP) (j £ (1,..., IV} \ {*})• The first 
order quantities being already dealt with in the last section, we need to bound those second 
order quantities in order to obtain bounds on (pi, 1 /*). For this, we compute their evolution 
equations along j th characteristics inside IZi ( j ■=/=■ i) using bi-characteristic coordinates 
(recall that yi is a parameter along Cj(zj)). 

First, notice that from (I22al) . we have for every j, k 6 (1,..., N} \ {*}, 

0 pj d Ai Afc pj pfs d \ k A i pj d 


(152) 


dyj A i A j dsi A i A, p k A* A^ Osi A j A i p k dy k 


We will use this identity to pass from (yi, ^-coordinates to (yi, 2 /*,)-coordinates. Also, 
(I22al I and (1451) imply 


(153) 


= Pi 

dyj A i — A j 




i)w k e k , 
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which we will use to compute the partial derivatives with respect to yj of functions of u. 
Finally, for any j ^ i, we compute using (122aD . (I521l . (HT1) . (lf>J)l> . 


(154) 


(155) 



V 


0 pi d 

- - V =- V' 

dy% Aj Aj dsj 


Pi d 
\j A i dsj 


(PjW J ) 


pipj (( 


Ei w j® i+ E 


7 V Q 


We are now ready to compute the evolution equations for and along a j th 
characteristic inside 7 Zi ( j ^ i). 

From 1331), GoJ). GD, we have 



Using then 1331), (IT531 m, we can write this as 


(156) 


d ( i) 1 

- Ti = -- 


dy, 


■ ( (El 4 wP ) (i) u^++ 4 (l 

\ V 


pi 


- 2 ' cL 


+ (a/-%(?p! (Ap A0( 


- (i) T p w p + ( V 


P^iiP^J 


_j_ (»W. 

+ (A,-A i ) 2 


where the inhomogeneous term 




^Mj, with 



















32 


D. Christodoulou and D. R. Perez 


depends only on first order quantities which are bounded (cf. (11151) . (11161) 1 inside IZi. This 
is immediate from (11551) . (1251) . (115411 . (11551) . for 


(157) 



qy£i p^i 


+ (^2( D ^Ajp)( X 9 - Xi)w P W q ^j Pi + (^(DeAji)^- 

P,Q Q 

p^i,qz^i q^i 

- (X c ’jp wP ) ((X) c A q ) pi + 

p \ q / 

P#* q^i 


- 4i ( (X 'Yp q w P 'w q ) Pi + (X^ 

\ P,q p 

P^i 

p^q 




where we have used EE}. ES}, ESI). 

Similarly, from (1551) . (11551) . (TUI) . and using (I54al) . (1551) . (1551) . we have 


_d_ 

dy. 


(i) j 
K J UJ J 


d ( 

• d 

dy,. V 

■ dyj 

9 1 

f P 

dyj \ 

i Xj 

C£) 

|co | 

1 

-(A, 


9 

(JL. 

dyj 

\dyt 


+ 


(Xj - Xi) 2 
PiP j 


(E 4 


7 V Q 
YpqW^W^ 


p,q 

p^q 


PiPj 


w p ) w j + ^ 7pqW P 

p,q 
p^q 


w p w J + 


Pi 


?(£*) 

p^j 

( X 4p wP ) wj + X 




p,q 

p^q 


7 J P q ) wPw<1 


7 P Q 
TpqW^W* 


P^iiP^J 


p,q 

p^i,q^i 

p^q 


)(&«)’ 


+ xfir, ((4”' + E 2 ]>’) 
'( E 


7 P q 

7 pqW^W^ 


p,q 

p^j^q^j 

p^q 


)( 


dy. 


+ 


Pi 


Xi- 


^((x^+x^^X 


7 V dyj 


P^J 


P^J 



P^i,Py^3 















On the formation of shocks of electromagnetic plane waves in non-linear crystals 


33 


With 051) . (1152ft . m, we write this as 


(158) 


d a 
dy■, 


(i) j 

y ’UJ J = 


( (E 7> P w") { \pi + (E 2 1ip wP ') 

\ p,q p 


( i ) i 
V J TjV 


P,Q 

v 

P9 £i, q9 £i 


P^i’P^j 


^ ((E(4+^H) (i W + (4+ 24 ) ( v^) 


P^i,P^3 


+ 


PiPj 


(\j — Ai) 2 

.2 , o-.J 


JP ^ JP (A p — Ai)(— W t p io p + W w p ))w j ' 


+ oTTfcF (£ T" , < A - - + <V)) 


(Aj Ai)" 5 Pp 

Py^Q 

p^i,q^i 
P^j^Q^J 


P 2 j 


(Ai - A0 : 


((? 


—^(A p — A i)(— (, V p w p + 
Pi k 


im p _L W^; p ) 1 v 




-I _^2_(Oiy* 

(Ai - Ai ) 2 


where 

















34 


D. Christodoulou and D. R. Perez 


Again, the inhomogeneous term ^ N J in (1158ft depends only on first order quant¬ 

ities which are bounded in TZi, since (1158ft . (1251) . (1154ft . (1155ft . (154a[) . (156[) . (OxTl) imply 
(159) 

= -(XI - dr)) ( f 

r 3 \\ x p P,q / ) 


P^iiP^j 


(E 


X r Xi r / j i 


)) ((^y+Z! 2 ^ 


p^i,q^£i 

Py£q 


i p w p )v l 


p^i 


+ ( (^ ( D e r Cj p ) (\r ~ A i)w P W r ) W 3 + ^ ( Ai r 7pq) (Ar - A i)w P W q wX pi 

\ p,r P,q,r / 

p^i, r z£i 

p^j r ^i,p^q 

+ f (^2 ( De Ai) ( A r “ A *) wr ) w3 +Y1 ( 2Z? er7ip) Or- - Ai)w P W r ^ 

\ r p,r / 


(55 ^p wP ) w °+ 55 'y 3 pi wP ' wq 


P^i,P^j 


p,q 

p^i,q^i 

p^q 


pz£i,r^i 

(j2drW )pi+CiiV‘ 


(djiW 0 + ^ 2 7 i r W r ) f (^2 7 vq wVwq ) Pi + (51 ( 2 ^P + dpW) vX . 


r^i pz£i,qz£i pz^i 

p^q 

We are now ready to show statement (Hull of Theorem [5] 


4.2. The bounds on the second order quantities. We first establish the boundedness 
of ^w' in R x [0,t*) \u fc 7 Zk- This, together with a bound on -§^pi in the same region, 
will yield bounded initial conditions at the boundary of IZi for ^V, and which are 

necessary to show the boundedness of the latter. This will then finally yield bounded 
inhomogeneous terms in the evolution equations along 1 th characteristics for pi and Vi 
inside IZi, making it possible to show bounds on those quantities, too. 

Our starting point is the evolution equation dun for -§^w l along an 1 th characteristic. 
For i G {1,... ,1V}, let Si denote the region in IZi where every j th characteristic (j ^ i) 
stays inside TZi for increasing times, i.e., 


(160) 5i = |(a:,i) 


G IZi 


Vj G {1,..., A} \ {i} we have: 


if Zj is such that Xj(zj,t) = x, then (Xj(zj,s),s) G IZi,Vs G [t,£*)|. 


Now define 


(161) TZ i =TZ l \S l . 

Then, for every point in TZi, there is a j G (1,..., AC} \ {i} such that the j th characteristic 
that passes through that point has to leave TZi before t*. As a consequence, we obtain 
from that 

(162) h"g) in«.. 

By m, dm and dS3- this implies that 

(163) w i = =0(1) in TZi. 

Moreover, by urn recall that w' = 0{iP) outside TZi. We conclude that, for all i G 
{I,---,AT}, 

(164) 


w' =0(1) 


in R x [0,1*) \ Si. 
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Now fix ( x , i)elx [0, t*) \ (J fc Sk = (R x [0, t*) \ (J fe IZk) IJ (U ; 'TZi'j , and assume i and 
Zi are such that Xi(zi,t) = x (i.e., Ci(zt) passes through ( x,t )). Then 


(Xi( Zi ,s),s) 


s e [o, t] | Pj ^|^J Sk'j — 0. 

k 


This means that we can integrate froii along Ci(zi) enjoying bounds on all the coefficients. 
Indeed, similarly to (I89al) . consider 


(165) 


V(t) 


max sup 

(x',t')g\^J k S k 

0 <t'<t 




Notice that V’ is non-decreasing (and, actually, even continuous) in t. By (11441) and the 
preceding argumentation, we have for all i € {1,..., N}, as well as for every x and almost 
every t such that (x, t) G R x [0, £*) \ (J fc 


d 

d i, 


d 

( d i, . \ 

dsi 

- W (x,t) 

< 

~dTi 

Vd~x W M 


< CV'(t) - I), 


where C and D are constants, depending only on a, S and the C^-norm of fo . We can 
thus argue entirely analogously to how we obtained m to conclude that 


(166) 


V\t)<e ct (v\ 0) + §) -§, 


where W(0) is bounded by the fact that fo is C K and has compact support in 7o = [—1,1]. 
Since t is finite by (fTTHi) . <nm we thus established that, indeed, is bounded for all 

i in the region R x [0, £*) \ [J fe Sfc. In particular, remains also bounded outside all 

the strips, as claimed in Theorem [5] (fTTIl) . 

Notice that, since pi = 0(1) by (11161) and (11341) . the foregoing observations and (11451) 
immediately entail the boundedness of -§^pi in R x [0,t*) \ (J fe Sk and thus in R x [0, t*) \ 
[J fe 7 Zk (Vi 6 (1,..., AC}). Note, however, that the bounds on and -§^pi get worse, 

the smaller we pick i9. Indeed, they blow up exponentially as •& goes to zero (see dm 

(firm and urn 

We now focus our attention on {^Vj }• and {}• inside 7 Zi. So fix i £ {1,..., AC}. 

If we want to integrate (11561) and HUD along a j th characteristic, we will need the initial 
values of and at the boundary of IZi (recall that yi serves as a parameter for Cj), 
i.e., along Ci(— 1) and C;( 1). But, by (I22al) . (TUll 


d _ d pj d 

dy-j ^ dx^X i — X j dsj ’ 


As a consequence, the bounds on -§^pj and w : ‘ just established, as well as (1521) . (1551) . 
(HMD , HMD and HJ imply by continuity that and are bounded on Ci(— 1) (J Ci( 1), 
for all j € {1,..., AC} \ {«}. For later reference, let us write 


(167) 


= max sup |^V 3 (a: , ,t , )|, 

} *' |JCi(i) 

^Ho = max sup |^\u 3 (a; , ,t , )|, 

(x'.tqec^-TlJCiCi) 


which are thus both finite. Similarly, define 

= max sup \^Tj(x' 

= max sup (x',t')\. 


(168) 
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Our aim is to show that both these quantities are finite. Fixing j £ {1,..., N} \ {*}, we 
have from (115611 and (11581) . respectively, that, along any j th characteristic inside IZi, 


A W T . = y ^Af\ p + y + {i) Mj 

OVi * ^ ^ 


(169) 


p^i 


_d_ a 

dy 


( V = y + y + {i 


:) A5- 


p^i pz^i 

where the coefficients are, by (11151) . (I116D and (11541) . at least 0($ 2 ), except for and 
^F 3 (p £ {1,..., AC} \ {i}), which are 0(1) and 0($), respectively. Set 

«A = 2 max sup I ^A p (x', t') I, 




(170) 


— 2 max sup \ i '^B 3 p {x',t')\, 

= 2 max sup |^M J -(a; , ,i , )|, 
■,-ji (*'>*')6'R i 

0® = 2 max sup \^E?(x',t')\, 

{xy F = 2 max sup I ^F 3 (a/, t') I, 

<*'■*'>€*« 

0f/V = 2max sup 1^AC J (a:', t') I, 

A. (i',t')6Ri 


which are all 0($ 2 ), except for and ^F, which are 0(1) and 0($), respectively. Then 
integrating along Cj yields 

(i) Tj < (l) G 0 + (iJ A (iy G + (iy B (iy H + {iy M, 

( V < {iy H 0 + (iy E (i) G + (iy F (iy H + (iy N, 


where we have used that, along cj n^> we have yi £ [—1,1]. Taking supremums over 
IZi, we obtain 


(i) G < (i) G 0 + (i) A (iJ G + (x W y H + (iJ M, 

(ijjj < (0# 0 + + (Qpdyjj + (i)]v. 


But since (i) A, (i )M, 0®, OW = 0(i? 2 ), while (iJ B = 0(1) and (iJ F = O(tf), and since (iJ G 0 
«ff 0 are bounded, we can re-write this as 

(i) G = 0( (i) Go + tf 2(i) G + (iy H), 

(iy H = 0( {iy Ho + tf 2(i) G + § (iy H), 
which, for sufficiently small •&, reads 

(i) G = 0( (i) Go + (ir ff), 

(iy H = O( (iJ H 0 + '& 2(iy G). 

Inserting the first line into the second yields 
(i) G = 0( (i ' ! Go + (iy H), 

{ir H = O( (iy H 0 + tf 2(i) Go + •& 2(iy H) = 0( (l) Go + 
for sufficiently small $. As a result, we obtain 

(171) (iy G = 0( {l) Go + (iy Ho) and (iy H = O( (iy G 0 + (iJ H 0 ), 

so that and ^uj 3 are bounded, as stated in Theorem El (lull). 
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At last, we are ready to show the boundedness of fii, v 1 inside IZi. In order to do this, 
we wish to integrate (fnni) and (fijii) along an I th characteristic inside IZi. We can re-write 
the evolution equations as 

d_ 

dsi 

d 


(172) 


~pi — A iii [Xi + Bf li V -{- , 


dsi 


v — A v i + B v i v + C „ 


with bounded (see ifira . (firm , (fim (irm i coefficients. Setting 

(173) Mi(t) = sup \iii(z'i,s'i)\ and N\t)= sup |i/(z-,s-)| 


(*'.«() 

0<s' <t 

which are non-decreasing in t, as well as 


z i eI 0 
0<s' <t 


(174) 


= sup 

\ A m (*'.*% 

(x' 


A v i = sup 

\A u i (x',t')\, 



B = sup 



B v i = sup 

\B„i (*',!') , 

(x' 

= sup 


(x f ,t')(E 7 Zi 


C v i = sup 

\C v t (x',t')\, 

(x f 



we obtain, for almost every t, 
d | | 

-\u l (zi,t)\ < 




ds 


dsi 




— ApiMi(t ) + B lli N (t) + C, 
< ~A v i Mi (t) + T3 v i N* (t) + ~C v i. 


Integrating along Ci and subsequently varying zt £ Io, we arrive at 


Mi(t) < A lli / M i (s)ds + B ll 


(175) 


N\t) < A v 


Mi{s)ds + B v 
a 


N\s)ds + C Ml t, 


N i (s)ds + C„it + N l ( 0), 


where we have used that pi(zi , 0) = ^ Pi{zi , 0) = 0, since pi(zi, 0) = 1. Note also that 
1V 1 (0) < V'(0) (cf. (11651) 1. To show that Mi, N l are bounded on [0, t*), we make use of 
the second statement in Lemma U for t < t* < To is finite (cf. (11181) . (11191) 1. 

This, therefore, establishes the rest of Theorem El (lulll. 


5. Application to electromagnetic plane waves in non-linear crystals 

This section is devoted to the application of the theory established so far to the case 
of plane waves in a non-linear crystal. We assume the material to be birefringent, non- 
dispersive, non-magnetic and exhibiting electric third order non-linearities in the energy 
density. We begin by recalling quickly the relation between the Lagrangian and the energy 
density in the electromagnetic theory, from which one can recover the Maxwell equations. 
After restricting to the aforementioned type of material and plane waves, we proceed to 
exposing under what assumptions our theory is applicable. Subsequently, we show the 
existence of Riemann invariants in a special case, showing thereafter their non-existence 
for generic third order non-linearities. As an experimental setup would require to deal 
with a material boundary, we end the section with a quick description of that situation. 
We then solve in full (albeit only implicitly) the special case above, and comment briefly 
on the problem under generic conditions. 
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5.1. The relevant equations. Inside matter, the Maxwell equations are 

divB = 0, divB = p, 


(176) 


4-S +curls = 0, 4-D - curl H = - J, 

at at 


where E , B, D, H, p and J denote the electric and the magnetic field, the electric and 
the magnetic displacement, the density of free charges and the density of free currents, 
respectively. These equations can be obtained from a variational formalism as follows (see, 
e.g., IChrOOI 'l. Let A = (^4q, Ai, A 2 , A 3 ) denote the electromagnetic vector potential such 


that, for i £ {1,2,3}, E x = and B l = y~) 


-Aie lls , where denotes 


the fully antisymmetric symbol (here, and in the following, indices are raised and lowered 
using the Minkowski metric diag(—1,1,1,1)). The field tensor F = dA, a two-form, is 

d_ a d_ 


given by F = 0 k F ^ dx>1 A ^ 


with Fpv = -£jtA„ - £ {0,1,2,3}, 


ar = t. We have 


/° 

-Ei 

-e 2 

-E 3 

Ei 

0 

b 3 

-b 2 

£2 

-b 3 

0 

Ei 

\e 3 

b 2 

-Si 

0 


= 


and we obtain the two equations on the left in (1 1 7611 from dF = d 2 A = 0. Given 
a (dual) Lagrangian density C* = C* {x,a,Q), where a and <1? denote one- and two- 
forms, respectively, we define the dual canonical momentum 'k* and the dual electro¬ 
magnetic displacement G*, both bivectors, by \k* = v _ Q |'l'* M1/ gf -p: A and G* = 
E^=o A respectively, where 

d 


>k* 


<9<E 


-C 


G* tl ' v (x) = 'S! i ’ tl ''(x,A(x),F(x)), p,v £ {0,1,2,3}, p ^ v. 


Setting 




and 


H' = -^-(C\-,A,F)), i £ {1,2,3}, 


we get 


(G* M ") = 


0 

D 1 

D 2 

- D 1 

0 

H 3 

- D 2 

—77 3 

0 

-D 3 

H 2 

—H 



The canonical momentum *k and the electromagnetic displacement G are then given 
1,T / M „cte M A dx u and G = i ._ 0 |G M „d:r M A dx v , respectively, with 


by * = 


antisymmetric symbol. We obtain 


(G^) = 


and G = 

,u=0 

/AIA - 

O 

II 

«< 

CO « 

w 

1 1 

2 '- T 

' 0 

H 1 

H 2 

—H 1 

0 

D 3 

-H 2 

-D 3 

0 

-H 3 

D 2 

-D- 


H 3 

-D 2 

D 1 

0 


Setting 


P(x)= (-q^£*)(x,A(x),F(x)), 
J\x)= (^---^-C k ^j{x,A{x),F{x)), 


i £ {1,2,3}, 


and defining the dual four-current I* by I*° = p and I* 1 = J 1 , i £ {1,2,3}, the current 

three-form 7 is given by I = 0 \l^dx x /\dx lx /\dx v , where I \= J] K =o 

The equations on the right in esj are then obtained from the Euler Lagrange equations 




(dG)* = /*, 


dG = I 
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which arise from varying A in the electromagnetic action <S[j4; T>\ = J C( ■ , A, F) of A in 
the space-time domain T> (cf. IChrOOl eq. (6.1.32)]). The current I being the canonical 
force l along a solution, I = (,(■, A, F), we write 


L 


3 

y ' -Lx^dx^ A dx 11 A dx” 


f'Xu 


£' 




* K 

i 



From the requirement that the action S'[.A; T>\ be invariant under gauge transformations 
A K» A + df compactly supported in the space-time domain T> , one shows (cf. IChrOOl 
p.265]) that l has to be a conserved null current , in the sense that dj = 0 and that 
+ g^-F v = 0 (see IChrOOl. §6.3&Ch.3]). Theorem 6.1 in IChrOOl then implies that 
l is, in general, of the form 

(o) (i) 

i = i(x, <£>) =C (x) + $(x)A C (*), 


where C, i e {0,1}, are closed exterior differential forms on space-time of degree 3 — 2 i. 
One concludes that the Lagrangian is of the form: 

(o) 

C(x, a, <E>) = £ (x, 'll?) + a(x) A i(x, $). 


The canonical electromagnetic stress tensor T* is defined through its components by 
(cf. IChrOOl §6.3]) 


T 


V 


( 0 ) 

(• ,F)S ^ 

A=0 


e (0,1,2,3}. 


Corresponding to a time translation field X, we have the energy-momentum density vec- 
torfield 

3 


p *! 1 _ 




u=0 


Taking X = we get 


Then P*° = £ is the field energy density and P* l = V 1 is the field energy flux (Poynting 
flux). From the above, we find in terms of E, B, D, H and slightly abusing notation that, 
along a solution, 


( 0 ) 

£(■ ,D,B) — E ■ D+ £* (■ ,E,B) and V{-, E, H) = E x H, 


where, now, 


E ~w e - 




( 0 ) 


If we assume £* not to depend explicitly on x° = t, the Maxwell equations imply 


= or", 


p ,=0 


fi—0 


or, 

%-£ + V-P = -E-J, 
at 

which expresses the differential energy law. 

We now focus our attention on electromagnetic plane waves inside a non-dispersive 
homogeneous dielectric material with no magnetic properties. Dealing with a perfect 
insulator, there is no free current, i.e., J = 0. For simplicity, we shall also assume the 
absence of extraneous charges, i.e., p = 0, so that we have 1 = 0. Note also that we assume 
to be well below the electrical breakdown of the material, i.e., we assume the fields to be 
weak enough so that the material actually behaves as an insulator. This is in accordance 
with our mathematical treatment, where we also showed that the smallness assumption 
on the initial data is preserved in evolution, even though infinite field gradients develop 
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in finite time. The Lagrangian density C* for a homogeneous material with only time- 
independent properties (hence our requirement of no dispersion) is invariant under both 
time and space translations. We can thus write, for some function L, 

C*(t, x, A, E, B) = L(E, B), 


and D, H are defined by 




h = M l • 


respectively. If the material has no magnetic properties, there is a function 5 such that 


L = 5(E) + -\B\ 2 


whence H = B and D = — i.e., D depends on E alone. Since the Maxwell equations 

cm?} are left unchanged by addition of a constant to L, we may assume that 5(0) = 0. 
In addition, we shall assume that —d5 is invertible around zero. In terms of the energy 
density (viewed as a function of ( D , B)) 


£{D,B) = E-D + L(E,B), 


we obtain, for i £ (1, 2,3}, 

■V=£f- 


<97> 




v-/ d 


;E’\D’+E'+^l 


-L 


dEi_ y 

=-Di 


dD 


-E j = E\ 


and 




d B' £ l^[dB‘ E:l ) DJ dEi L ) dB' EJ + dB* L H ’ 


as expected. Since D depends only on E, we can write 

£(D,B)m<5(D) + ±\B\ 2 , 

where 

<&{D) = E-D + ${E). 

Note that 5(0) = 0 entails ©(0) = 0. 

The aim of what follows is to study plane wave solutions of the Maxwell equations for © 
given by its third order Taylor expansion around zero. For emphasis, we shall henceforth 
write x = x 1 , y = x 2 , z = x 3 , and, correspondingly, 

D x = D\ Dy = D 2 , D z = D 3 , E x = E\ E y = E 2 , E z = E 3 , 

B x = B 1 , B y = B 2 , B z = B 3 . 


An electromagnetic plane wave travelling in the rr-direction is characterised by the fact 
that the field variables depend only on x and t. From the constraint equations divB = 0 
and divZ) = 0, we obtain that D x and B x are constant in space. Their constancy in 
time is a consequence of the two other equations in inm We may therefore assume 
that D x and B x vanish throughout. The remaining equations governing the evolution of 
{D y , D z , By, B z ) reduce to 


D 


(177) 


4- 9 R 

y+ di Bz 


d_ 

dt 


d_ R ^d_ F 

dt By dx Ez 


= o, 
= o, 


D z - 


d_ 

dt 




i s -= 0 ’ 

+ Ae» = 0, 


where E = ^ In the type of material under consideration, modulo the correct 

alignment of the rectangular coordinate axes, we can assume that, up to terms of fourth 
order, 

©(D) = \k!DI + \k 2 D 2 z + CmDl + C 112 D 2 y D z + C 122 D V D 2 Z + C 222 D 3 z , 

where Ki,K 2 , C\n,C\i 2 , C\ 22 ,C 222 £ R are constants (note that there is no need to con¬ 
sider terms linear in the components of D, as they would not alter the Maxwell equations 
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(11771) ). In order to obtain a strictly hyperbolic system, we restrict ourselves to birefringent 
materials, i.e., media where the refractive index depends on the polarisation of the light. 
This situation occurs in certain crystals and translates into the condition 

(178) 0 < AT < AT < 1. 

The energy density is 

(179) £ = \B 2 y + ^B 2 z + ^K 1 D 2 y + ^K 2 D 2 z +C 111 D 3 y + C 11 2D 2 y D z + C 1 22DyD 2 z + C222D 3 z . 

We then obtain from E = ^ ^, 


(180) 


E y = K\D y + 3CmDy + 2C 112 D y D z + C\ 22 D 2 Z , 
E z = K 2 D z + C U2 D 2 y + 2C 122 D y D z + 3C 222 D 2 . 


We now show that (1TT71) . taking into account cm falls in the framework of our theory. 
More precisely, we prove 

Lemma 5. Let u = (D y , D z , B y , B z j. Then there is a matrix function a = a(u ) such 
that (HTHl for E given by ifTHUl) is equivalent to 


(181) 


d , J 
-u + a(u)-u = 0 . 


Moreover, under the assumption mm there is a 5 > 0 such that the system cm is 
uniformly strictly hyperbolic on B$( 0). Finally, (11811) is genuinely non-linear, whenever 
Ci n and C 222 are non-zero. 

5.2. Proof of Lemma[5j Let u = (D y , D z , B y , B Z ^) T . Then the first row of (11771) is 

9 1 d 4 d 2 9 3 

di u + ^ u = 0 ’ di u =a 

Using GBB , the second row of mu translates to 

£«» - K 2 ^-U 2 - 2Cll2U 1 ^-U 1 - 2C'i22U 1 ^-u 2 - 2 Ci22« 2 ^-u 1 - 6C 222 u 2 ^-u 2 = 0, 
ot ox ox ox ox ox 

and 

9 4 , ry 9 1 ayi 1 9 1 1 9 2 | r\/—i 2 9 1 r\/~i 2 9 2 „ 

—u + Ki—u + 66111 U —~u + 26112 U —u + 26112 U —u + 26122 U —u =0. 
ot ox ox ox ox ox 

Thus, setting 
(182a) 

(182b) 
and 
(182c) 

we see that the Maxwell equations mu are equivalent to 


d\ = di(u) = Ki + 66 'inu 1 + 2 C 112 U. 2 , 
d 2 = d 2 (u) = K 2 + 26'i22U 1 + 6C 222 u 2 , 


C = c(li) — 2C112U 1 + 26l22« 2 , 


9 d 

-—u + a—u = 0 , 
ot ox 


with 


0 0 

(183) 0 *“(“)= |-c(«) -&(«) 

v di(u) c(u) 

We have thus proved the first statement in Lemma [U 
Now set 

di + d 2 


0 IN 
-1 0 
0 0 
0 oy 


(184) 


m = m(u) = 


and 


r = r{u ) = 


di — d 2 


so that m + r = di and m — r = d 2 . Notice that 
(185) di( 0 ) = AT, d 2 {0) = K 2 


and c( 0 ) = 0 , 
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whence, setting mo = m(0), rg = r(0 ), 

,, Ki + ^2 

(186) m 0 = --- 


and rg = 


A i - A' 2 


Pick then any h > 0 such that 


*2 




(187) 
so that 


h < 


min|r 0 , mo"—ro , 1 — (m^+ro) | 


2mo 


ro 

AT 1 

2mo 

2m 0 J 

rg 

i + a 2 

2mg ’ 

2mg 


Then it is easy to check that, by (fim 

0 < (1 - h)K 2 < (1 + h)K 2 < (1 - h)K 1 < (1 + h)Ki < 1, 
for 

a+ h) Kl < Ki S x + 5 a) = i - Ka S 7 5° < i. 


A i + A 2 


AT + AT 


(1 - h)K x - (1 + h)K 2 > AT - A 2 - A " ir ° + Azr ° = 2r 0 - r 0 = r 0 > 0. 

2mg 

Let now <5 > 0 be so small that 

|AT - rfi(v)I < hKi, |AT - d 2 (v)| < AA 2 , 


and 


c(i>) 2 < min{ Ai (A' 2 - h(Ai + A 2 )) , (1 - A 2 )(l - AT - h(Ki + A 2 ))j, 


(188) 


for all v € -B 2 i(0) = {?/ 6 R 4 | |t/| < 2<i}, and let u G B^si 0)- Then, by (118411 . (118611 . 
di + d 2 AT + A 2 di — Ai d 2 — K 2 


+ ■ 


+ 


/AT+AT fcAi hK 2 I<! + K 2 hKx hl< 2 \ , , 

£ {— - 2 -+ — + — ) =((l-*W(l + ft)mo), 

so that 

m > (1 — h)mo and 1 — m>l — (1 + h)mg. 

Since (118711 implies h < , we have by (11781 , 

. rg Ki+ K 2 I<i — K 2 1 „ 3 

(1 - h)mg > mg - — = --- - - = -Ai + -K 2 > 0 


and 


,i , „ , ro Ai + A 2 AT — A 2 3 1 

(1 + h)mg < mg + — = ---+--- = -Ai + -A 2 < 1, 


As a result, 

Define 

(189) 


(1 — h)mg > 0 and 1 — (1 + h)mg > 0. 

0 < m < 1. 


R = R(u) = \J r 2 + c 2 , 
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and observe that, by (11841) . (11881) . (11861) . 

R 2 = r 2 + c 2 < 2 + Kl (K 2 - h(Ki + K 2 )) 

= Tq + 2 hr 0 m 0 + h 2 ml + (mo + r 0 )(m 0 — r 0 — 2/imo) 

= ((1 — h)mo) 2 < m 2 . 

Similarly, (11841) . (I188D . (11861) yield 

R 2 = r 2 + c 2 < 2 + (1 - K 2 )( 1 - Ki - h(Ki + K 2 )) 

= rl + 2hr 0 mo + h 2 ml + (1 — m 0 + r 0 )(l — m 0 — r 0 — 2hmo) 

= (l — (1 + h)mo) 2 < (1 — m) 2 . 

Finally, (118411 . (11881) . (118611 tell us 

„ m, o - . Ki-k 2 a'i + k 2 

R = r 1 +c 2 > r > - - - h - - - =r 0 — hm 0 , 

whence, by (fl87l) . 

m — R < (1 + h)^— — R < mo — ro + 2hmo < mo 


and 


m + R> (1 — h) K 2 + R > mo + ro — 2/imo > mo. 


We have thus showed that 


(190) 0 < m — R < mo < m + R < 1. 

We are then ready to compute the eigenvalues of a. Writing (cf. (11831) . (11841) 1 

/ 0 0 0 i\ 

_ | 0 0 -1 0 I 

a —c —m + r 0 0 I ’ 

\m + r c 0 0/ 

we calculate 

det(A id —a) = A 4 — A 2 (m — r) — A 2 (m + r) + (m 2 — r 2 — (?) 

= (A 2 ) 2 - 2mA 2 + (m 2 - R 2 ) 

= (A 2 ) 2 — ((m. + R) + (m — R))\ 2 + (m + R)(m — R) 
= (A 2 — (m + 7?)) (A 2 — (m — 7?)), 
and obtain the eigenvalues A of a: 

(191) A = A(tt) = ±Vm ± 77. 


Notice that these are, by (1TM1) . real and of modulus strictly smaller than one. Moreover, 
the four values are distinct, uniformly on Bf(0). Indeed, by (11901) . we have on 77^(0), 


(192) — 1 < —\Jm + R < -/mo < —y/m — R < 0 < %/m — 7? < /mo < %/m + 77 < 1, 


which, therefore, holds uniformly on T?f (0). As a result, the second statement in LemmaO 
is proved. For the rest of this subsection, we shall assume that u £ 7?|(0). Notice that, 
for u = 0, we obtain by (11841) . (11891) . (11851) . 


(193) 


A(0) 


±y/K~ 2 


Also, it is easy to see that (see 


GED and (11851) 1 



( 1 \ 


f o \ 


f o \ 


0 


1 


1 

eoi = 

o 

, eo2 = 

-s[K 2 

, e 0 3 = 

s[K 2 


wmj 


\ 0 ) 


\ 0 J 


and eo4 



(194) eoi 
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are eigenvectors of a(0) corresponding to the eigenvalues Ai(0) = \JK i, A2(0) = VK 2 , 
A 3 (0) = —\TK 2 and A4(0) = —yJK 1 , respectively. Now note that the basis {eoi}ie{i,..., 4 } 
is orthonormal with respect to the scalar product 


(195) 


(v,w)o = y^gijVW 3 , 


V, W £ 


where gij, i,j £ {1,... ,4}, are the components of the matrix 

(196) ff = diag (^’2^’2^)' 

We will use the norms || ■ ||o and || • ||o* induced by this scalar product on R 4 and R 4 , 
respectively, to normalise the eigenvectors of a(u) and measure their duals (cf. Subsec¬ 
tion 12.11) . 

For c ^ 0, let 

(197) g, = — -— = — (r ± 7?) and ju = — —-— = - (r =F R). 

We have 

(198) 

One easily checks that 


/r/r = 4-(r 2 — 7? 2 ) = — 1 and /r + /Lt=—. 

c 2 c 


(199a) 


^A — & n 



where 

(199b) = sign(M)5 A>±v ^+fl + S x,±V^=R e { — 1,1}, 

is an eigenvector of a for the eigenvalue A with dual 
(199c) e* x = 

Having 

ll e *llo = {e\,e\)o = g ^1 + g 2 + ^ ) 


and 


H'~VrA 11 2 / ★ A -~*A\ 

11^ No* )o* 


2 + (2A ? + t^ + ^Y 


k 2 


8A 2 (1 + g 2 ) 2 

^(A 2 + AV + lf2 + W), 


A'i 


we normalise 

( 200 ) 

so that 

(201) UcaIIo = 1, 

and 


e\ = 


1+jJ 2 I X 2 ( _1_ 1 \ 

2 1" 2 V K 2 ^ K 1 I y 



+ ^ (k 2 + KT) , ^ 

2A(l + /x 2 ) 1,M )’ 


* A c A 

e ey = 0 \r 


e 0 * = 


/( 


1 + M 2 + A 2 ( 


+ M 2 ) + A 2 + Aip 2 ) 


2|A|(1 + fi 2 ) 
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for 


0><V( A + A')(l + fin') 


ev= 2A(l+V) ( Aw+A + A + A ^) =- 2A(1 + fj, 2 ) 


and, in case A ^ A', we have either /r = \x , whence A = — A', or fi = //, whence fifi' = — 1 
by (11981) . 

On the other hand, for c = 0 (which occurs, e.g., when u = 0, or when C 112 = C 122 = 0; 
cf. fl 1 82cfl 1. a becomes 

/0 0 0 

0 0-10 
0 -d 2 0 0 

\dt 0 0 0> 

and, since R = r in this case, we obtain 


A = 


± \fd[ 
± \fd~2 


Now the problem decouples into a pair of two-dimensional equations for (u 1 , u 4 ) and 
(u 2 ,u 3 ), respectively (i.e., for ( D y ,B z ) and ( D z ,B y ), respectively). The (normalised) 
eigenvectors of a are easily seen to be given by 


1 




e ±v^ 


( 1 
0 
0 


*±y/d\ _ 


1 d l 

2 T 2Kl 


1 d l 

2 ' 2Ki 

i2\/ di 


1 _ 1 _ ^2 

2 ^ 2 K 2 


( 0 

+V<h | ’ 
V 0 


'I _L <L 2 

±2-\/d2 


(±\/di, 0, 0, l), 


(0,±v/dl,-l,0). 


We show that our eigenvectors mm for 0 can be extended by continuity to c = 0 
and coincide there with the appropriate eigenvector from above. Let 

(202) /x± = - (r ± 7?) = - fr ± \J r 2 + c 2 

For 0 < |c| < r, we have 

r -¥^M(T)y 


Thus 

so that 
(203) 


M± 


2r /1 cl \ 

fl+ = ~ + °\v) 


and 




&)' 


and sign(/r±) = sign(c) for |c| small enough. A short computation then gives 


++ 


\A + M+ 


li- 


= sign(M+) + C»(y), 


\A+ 

As a consequence, 


= O 


(*)• 


\AT74 

1 

\/l + M- 




= O 

= 1 + 0 


(£)• 


c—>-0 

^±\/m+R ^ ^ 

c—>0 

e ±Vm-fl ^ e 


±V ’ 


giiV'm + fl e ~ > ( 1 e *±\f+ 
*±Vm-R e~>9 *±\fd2 


which proves the claim. 
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Let us now compute the “gradient” of the matrix a. Let u = (ii 1 , u 2 , w 3 , m 4 ) T be a 
variation of u. Since the coefficients of a are affine functions of the components of u, we 
obtain simply 


( 204 ) 


where 

( 205 a) 

( 205 b) 

( 205 c) 

It follows 


( 0 0 0 0\ 

0 0 0 0 | w 

-d 2 0 O ’ Vs G 

d\ 600 / 

d\ = 6C111M 1 + 2C112U 2 = 2 (3C111U 1 + C11211 2 ), 
d2 = 2C122M 1 + 6C222U 2 = 2 (Cl 22 M 1 + 3 C I 222 'U 2 ), 

6 = 2C112M 1 + 2Cl22'6 2 = 2 (C'll 2 U 1 + Cl22U 2 )- 


*A / r-, \ 


l + M 2 + A 2 (^ + £) 


0 

(A M ,A,-l,/r) I _ , °_ 


+ m ' 2 + A' 2 (^ + ^) 

1 + /r 2 + A 2 + j^) pp'dr + ri 2 + (p + //)c 

\| 1 + M' 2 +A' 2 (^ + ^) 2 A (1 + /r 2 ) 


fj,' c — d2 

fi'di + c 


As a result, we obtain for u = e A ", 


CVA" = CVA"(«) = e * X ( u ) (u) a ( U )) e A' («) 


= <W<V, 

/r/i'( 3 Cin 

o 1 + M 2 + A 2 (-^ + f r ) 

' (1 + ^ + A*(£ + £)) (l + ^" 2 + A" 2 (^ + O) 

n" + C112) + (C\2ili' + 3C222) + (/r + /Lt , )(C'll2// , + C122) 


A (1 + At 2 ) 

( 206 ) = <T M oyoy/ 

o l + M 2 + A 2 (-^ + -| r ) 

' (1 + ^ + A' 2 (i + £)) (l + M" 2 + A" 2 (^ + <£)) 


3 C l iii/r/i , /r” + C , ll2(/A/r , + mm” + I-”) + Ci22(/r + + /r”) + 3C222 

A(l + P 2 ) 


Notice that c A /, A / = c x , x „. Notice also that at least two of the /r, //,//' are identical. I11 
particular, we have 

3 <t 3 

( 207 ) e AA =- . = (Cm/i 3 + C112M 2 + C122M + C222) ■ 

A(1 

As a consequence, we have for 

A = ±v^hTR : cL f ;Cm, 

±V^V 2 + 2kT 


and for 

A^iV^r^R : f C222, 

which shows that the system + a-^u = 0 is genuinely non-linear around u = 0, as 
soon as Cin,C222 7^ 0 . This establishes the last statement in Lemma [ 5 ] completing the 
proof. 
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Remark. If we want to ensure that c AA < 0 for all sufficiently small |u|, we replace e* x 
and ex by o\e* x and a\e\, respectively, where 

( 208 ) <J\ — sign(Cin) sign(A)< 5 A>±v ^ T7 j - sign(C' 2 2 2 ) sign^^^^rR e {-1,1}. 
(Note that, by definition, o-\ = —cr\.) We then have 


( 209 ) 

A 

C A'A" 


0 > 0 > ,C ’V" a >^ a y a \" 

O 

1 +p 2 + A 2 (^- + 

A(l + p 2 ) \ 


(l + p' 2 + A' 2 (^ + £)) 

(l +M "2 + A"2(-l- + *£)) 


■ ( 3 Cu_ipp /r 4 - Cn 2 (/r/r + /r/r T- p /r ) + Ci22(p + M T h ) “t - 3C222) 


and 


( 210 ) 


A 

CAA = 


o _3 _3 
^ a u a X 


*(i + » 2 )]/W+f(ik + W) 

Remark. Observe that, for u = 0, we have 


(Cm/x 3 + Cll2/i 2 + Ci 22M + C222) ■ 


( 211 ) 


min c AA (0) = min< 


max c AA (0) = maxi 


f 3 

3 , ,1 

{^rr |Cnl1 

’ VKl 102221 '/ 

f 3 , 

3 , ,1 

{vxr 1 11,1 

5 xm 102221 j 

max A |c£ A (0)| 

if and only if 


|Cm| 

IC222I 

which, by cm, is strictly larger than one. 



This finishes the discussion that our theory is applicable to plane electromagnetic waves 
in a birefringent crystal with no magnetic properties and cubic non-linearity in the energy 
density, assuming no dispersion. Except for a smallness condition on the field we note 
that we only need to require Cm, C222 7^ 0 . We proceed by observing in the next two sub¬ 
sections that there exist Riemann invariants in the decoupled case c = 0 , thus simplifying 
the problem considerably, whereas no such simplification occurs in the case c ^ 0. 


5 . 3 . On the existence of Riemann invariants in the case c = 0 . Let us assume that 
e = 0 (i.e., we assume that C112 = C122 = 0 ; cf. (U 82 cD h We already observed that this 
implies the decoupling of the two polarisations. Now we show that this actually implies 
the existence of Riemann invariants. 

First, notice that C112 = C122 = 0 implies that di(u) = difu 1 ) and d 2 (u) = d2{u 2 ) 
(cf. IU 82 al) and (I 182 bl) 'l. So, indeed, the two decoupled systems do not depend on one 
another. We can then define 


( 212 ) 



It is easy to see that a becomes diagonal in these coordinates. Indeed, we have 


/ x/d4yPj 0 ° l\ 

0 x/d 2 (u 2 ) -1 0 

0_ - x /d 2 {u 2 ) -1 0 

\— \/di(w 1 ) 0 0 1/ 


dm(u) 
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so that 


(dm(u)) 1 = i 


( y/di^u 1 ) 

0 

0 

V 1 


0 


y/d -2 (u 2 ) y/d 2 (u 2 ) 

-1 -1 

0 0 


y/ diO 1 )^ 

0 

0 
1 


and therefore 


dm(u)a(u) (dm(u)} 1 = 


^ di(n 1 ) 0 0 

0 yj d 2 (u 2 ) 0 

0 0 —y/d^vF) 

VO 0 0 -y/MvPjJ 


0 \ 
0 

0 


We thus obtain that + a(u)-j^u = 0 is equivalent to the system 


d i / , , n 5 i 

+v / dl(w i)_ m =0, 

d 


^.m 2 + y/d^)gm 2 = 0, 


= 0, 


a 


-m 4 -y^)-m 4 = 0, 


which we write more compactly as 

9 


(213) ^m i(1 - j)+{5 - <)j ' + (-l)V*(« i )^"» i(1 - J ' )+(5_i)i = 0, 


>6{1-2}J6{0, I}- 


We still need to express u 1 and u 2 in terms of m 1 , m 2 , m 3 and m 4 . From (I182al) . 
(1182bf) we have 

16(1,2}, 

whence 


(214) / y/di(s)ds = 


yfdFivF) = y/Ki + dCiuv?, 

(Ki + dCuiU^ — (Ki)i 


9 Ci, 


\fKiu\ 


, if Cm y 0, 


if Cm = 0, 


i6{l,2}. 


But from (Em we see that this also equals -(ml — m 5 l ). As a result, we obtain for 
16(1,2} 

(215) ii = u 1 (m l , m 5 ” 1 ) = 


(Ffi) 2 + ^(m 1 - m 5 


— Ki 


6Cm, 

:(m‘ — m 5_l ), 


if Cjii ^ 0, 
if Ciii = 0. 


2yfKl 

Inserting this back into (I182al) . (1 182bl) and then into (12131) . we obtain for i £ {1,2}, j £ 

(0,1} 


d m 

- f~ 

dt 


2ij+5j 




d 

dx 


i-2ij+5j _ 


0, if Cm y 0, 


|m’- 2,i+5j + (-1 ) j VlC-^m i - 2ij+5j = 0, 


if Cm = 0. 

The first equation being compatible with the second, we have for every Cm £ R, 


(216) wr 


& _ i — 2ij+5j 


+ (-iy[(K^+^l(m i -m 5 - i )) ^ 


i-2ij+5j _ 


= 0 , 

*e{i,2},je{o,i}. 


We conclude that m' 2lJ+5 -' is constant along the (i — 2ij + 5j) th characteristic with speed 


(217) 


S 9 C, 


Ai-Mi+Bi = (-1) J ( {KiY* + - m 5-i ) 
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As a result, the characteristics must be straight lines as soon as we are outside a region 
where characteristic strips overlap. 

Next, we want to investigate whether or not a similar argument can be performed in 
the case c ^ 0. It turns out that the answer is negative. 


5.4. On the non-existence of Riemann invariants in the case c ^ 0. Notice that, 
in general, the existence of Riemann invariants for an IV-dimensional (N £ N \ {0}) 
quasi-linear system 


d m 3 

—w + a(u) —u = 0 

depends on the existence of coordinates m = m(u ) in R lv such that the rows of dm are 
(each separately) proportional to the elements of a basis {e* x }fL 1 of eigen-l-forms of a 
(cf. Remark [I] on page [4]), i.e., there must exist functions /u,i such that 


[iidm l = e ", i £ {1,... ,N}. 


We have the necessary integrability condition 

de* 1 A e** = (dyt A dm 1 + p id 2 m 1 ) A (/ udm *) = fiid^i A (dm 1 A dm 1 ) = 0. 

In order to check this for the case at hand (i.e., N = 4, a given by (11551) and c ^ 0), 
we compute 

de* A = (d u i A — 9 u 2(A/r))dM 1 A du 2 + {—d u 3 (Xf-t^du 1 A du 3 + (d u i fi — d u i(\fi^du 1 A du 4 
+ (—d u 3X)du 2 A du 3 + (d u 2 fi — d u A\)du 2 A du 4 + ( d u 3fi)du 3 A du 4 , 


where we take 

e* A = Xfidu 1 + A du 2 — du 3 + /xdu 4 

for the dual eigenvectors of a (compare with (1199c!) '). Since, by (119711 . (11911) . (11891) . (11841) . 
(I182al) . (I182bl) . (I182cl) . the quantities A and fi are independent of u 3 and u 4 , this reduces 
to 

de* A = (<9 u iA — d u 2(Xfi))du 1 A du 2 + (d u if^du 1 A du 4 + (d u 2fi)du 2 A du 4 . 
Therefore, 

de* A A e* A = (— d u i A + d u 2 (A fij'jdu 1 A du 2 A du 3 

+ i A — (Ap.) — A<9„i/r + X^,d u 2 /r) du 1 A du 2 A du 4 

+ ( 9 u i fi)du 3 A du 3 A du 4 + ( d u 2fi)du 2 A du 3 A du 4 , 

which has to vanish. Comparing coefficients, we obtain that fi needs to be constant and 
that A has to fulfil 


(218) d u i\ = fid u 2\. 

Since, by the first equation in (ITMll . fi being independent of u implies that is independent 
of u, we obtain from the second equation in (1T551) that ^ has to be a constant as well. On 
the other hand, considering points closer and closer to the trivial state, (12031) implies that 
/x_ = 0 and = o. Hence, again by (15551) . f = 0, implying that c = 0. But this is a 
contradiction to our assumption that c ^ 0. As a result, no Riemann invariants exist in 
the case c ^ 0. 

In the remaining part of this section, we want to comment on a more realistic setting 
one might encounter in experiments. 


5.5. A sketch of an actual experimental setup. If one were to set up an experiment 
for comparison with the theoretical study conducted thus far, one would need to take into 
account at least a material boundary. In what follows, we shall assume that our crystal 
extends infinitely in the half-space |(*,j/,z) £ R 3 | x > 0} and that a suitable plane wave 
source is located in the plane {x = a?o} for some xo < 0. For simplicity, we shall assume 
that the half-space {(a;, y, z) £ R 3 | x < 0} is vacuum, although it should be easy to adapt 
what follows to the case of, say, air. The initial conditions for the electromagnetic field are 
obtained from a smooth light profile emanating from xo and compactly supported in time. 
We will later choose the time origin so that the light first hits the crystal at ( x,t ) = (0, 0). 

In order to deal with the material boundary, we obviously have to set up jump condi¬ 
tions at x = 0. Those are, as usual, obtained from integrating the Maxwell equations over 
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suitable regions (infinitesimally flat pill boxes, respectively loops, based on arbitrary sub- 
domains, respectively curves, on the boundary) and using Gauss’ and Stokes’ theorems. 
In general, one obtains for insulators (at rest) without extraneous charges and currents 
that the normal components of B and D, as well as the tangential components of E and 
H, have to be continuous across the boundary. We denote this as 


(219) 


[B ± ] = 0, [D ± ] = 0, 

[£,|]=0, [ff|,] = 0. 


For plane waves travelling in the ^-direction, the boundary being {x = 0} here, the top 
line has no effect other than ensuring that if we assumed B x = 0, D x = 0 initially, the same 
holds on the other side of the boundary. The bottom line, on the other hand, translates 
into 


(220) [E y }= 0, [E z }= 0, [tf w ]=0, [H z } = 0. 

These conditions, of course, have to be converted into jump conditions for D and B, using 
H = B on both sides, as well as E = D if x < 0 and (fTHUl) if x > 0. We observe that 

d(Ey,E z ) fd! c \ 

d{D y ,D z ) \c d 2 J ’ 

so that ( D y ,D z ) can indeed be expressed in terms of (E y ,E z ) in a neighbourhood of the 
solution (D y ,D z ) | (J5biJ5i)=(0i0) = (0,0). 

We now give the complete solution of this situation in the decoupled case c = 0, i.e., 
when Cn 2 = C 122 = 0. 


5.6. The decoupled case c = 0. Both inside and outside the crystal, we take ifTTSl) as 
an ansatz for the energy density. In the region outside the crystal, i.e., when x < 0, the 
constants Cm, G112, G122 and C222 all vanish, and K 1 = K 2 . For simplicity, we shall 
assume that we are in vacuum, so that K\ = K 2 = 1 . Inside the crystal, we assume 
<nm C112 = C122 = 0 and, for simplicity, we only consider the genuinely non-linear case 
Cm ^ 0, G222 ^ 0. Let u = (D y , D z , B y , B z ) be a solution of J|u + a-^u = 0, where 
a is given by GS 3 ) with the aforementioned values of the constants in accordance as to 
whether x is positive or negative, and such that the jump conditions (12201) are fulfilled. 
In what follows, we describe this solution in terms of Riemann invariants. We can use 
Subsection 15.31 to do this in both the vacuum and the crystal regions. To avoid additional 
sub- or superscripts, we shall denote the Riemann invariants inside the crystal by m, but 
those outside the crystal by n. Let us start with the vacuum region (x < 0). Here, 
since the equations are linear, we expect the solution u to be the sum of an incident 
wave up) travelling in the positive ^-direction, and a reflected wave up) travelling in the 
negative ^-direction, so that u = up) + up). Given an initial two-dimensional profile 
f € Gf (R;B|(0)) for (D y ,D z ) = (u\u 2 ) with compact support in ( xo — 1,2:0), we can 
parametrise the incident wave up) by 

/ f 1 {x 0 + x - t) \ 

(221) wp)(apt) = (*o,0), 

V f 1 (xo + X - t) ) 


as can be easily checked using (1X771) . Here, 5 is chosen small enough according to Subsec¬ 
tion for the crystal region. Notice that we chose the time origin in such a way that 
t 1 —^ lim uii)(x,t) is compactly supported in (0,1). The Riemann invariants np) related 

cc/*0 

to this wave are (cf. ( 12121) 1 

up) T up) \ 

^p) U(p I 
Ufi)^ arp) I 

^(i) T tap) J 



/2f (xo + x — t)N 
2 f 2 (a:o + x - t) 

0 
0 


( 222 ) 


np) (x,t) 


x G ( 10 , 0 ), 
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which is consistent with uu) travelling in the positive ^-direction. Similarly, we can make 
the following ansatz for the Riemann invariants ri( r \ representing the reflected wave u ( r ), 


(223) 


Tl(r) [x,t) 


s 1 

2 S 2 (x,t) ’ 

\ 2 S 1 (x,t)J 


x < 0 , 


for two unknown functions g 1 and g 2 (notice that, by linearity, n = nu\ +n( r )). Since the 
latter are constant along the characteristics with speed — 1 , they must depend solely on 
x + t. Slightly absuing notation, we then get for the reflected wave ut r \ using (12121) . (I215D . 

/-g\x + t)\ 

(224) u (r) (x,t)= , x<0, 

\ g 1 (x + t) ) 


which is consistent with the Maxwell equations CEB- We shall promptly see that the two 
unknown functions g 1 and g 2 are obtained from the jump conditions (Em We are going 
to need the following expressions for E y , E z . H y , H z in terms of n 1 , n 2 , n 3 and n 4 , which 
are readily obtained from (12121) . (12151) . using that, in vacuum, E = D and H = B. We 
have, for x < 0, 

E y = Dy = u 1 = Un 1 - n 4 ), E z =D z =u 2 = i(n 2 - n 3 ), 

(225) - f 

Hy = B v = u 3 = —-(n 2 + n 3 ), H z = B z =u 4 = -(n 1 + n 4 ). 


Regarding the region inside the crystal, we consider the Riemann invariants m = 
m(x,t) (x > 0) representing the transmitted wave. Since m 3 and m 4 are constant along 
characteristics with negative speed (see (BUD ), and assuming that the incident wave from 
above is the only source, it is clear that they must vanish, as every characteristic C 3 or C 4 
intersects points in the trivial region. We thus have 


m(x,t) 


i 1 (x,t)'< 
m 2 (x,t) 
0 
0 


x > 0, 


and observe that, by ET51) . (Em . the characteristics in the crystal region are all straight 
lines. We are going to use Em to obtain m 1 and m 2 . Before doing so, however, we need 
to express m 1 , m 2 , m 3 and m 4 in terms of E y , E z , H y and H z , inside the crystal region 
x > 0. We have that 


U 1 = Dy 


y /K\ + 12Cing„ - K 1 
6C111 


and 


u =D Z 


y /Kj + 12 C 2 2 2 E z - Ka 
6C222 


are the unique solutions of Em that vanish when E y , respectively E z , vanish. We then 
obtain from Em . Em . and using H = B, that, for x > 0, 


(226) 


r 4 , (RTi+ 667111^)1 -K? „ , {Kl + 12Cm^)i - k\ 

m =u + - ochTi - = Hz + -9^r- 1 3 

2 3 (K 2 + 6 C 222 U 2 )i - kI (Kl + 12C222E z )i- k! 

m = ~ U + - 9 ^- = - Hy + - 9C222 - 

3 3 (K 2 + 6 C 2 22 U 2 )§ - kI {K% + l 2 C 222 E z )i- kI 

m — u — ft ij 

9C222 9(^222 

4 4 (#1 + GCiiiu 1 )^ - k\ {Kl + 12CmEy)i - k\ 

m =U -9^1- =Hz -9^1-• 


We are now ready to use Em to solve for m 1 and m 2 , as well as for n 3 = n ( r ) 3 = 2 g 2 
and n 4 = n ( r ) 4 = 2 g 4 from above. 
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For * 6 {1,..., 4}, let 


no(t) = lim n(x,t ) = 
xyo 


^fixo -t)> 

2f(*o-t) 

20 2 W 
20^) 


mo(t) = lim m(x,t) = 
x^O 


j lim m 1 (a;,t)\ 
lim m 2 ( x, t ) 

:c\*0 

o 

o 


Using 


and dm we obtain 


1 

n 0 + n 0 1 

m 0 = 

—-u 

2 

2 

no + n,Q 

m 0 = 

-b 

2 

3 

n o + n o 

m 0 — 

2 

4 

n 0 + n 0 

to 0 - 

2 


9Cin 

2 3 . 3. 


9C*222 

^>2_ 3 , 3 


9C*222 

(A? + 12Cin * - k\ 


Using 


m o(t) = + 

ml{t) = J 2 (t) + 0 2 (t) + 
0 =f(t) + s 2 (t)- 
0=f(t)+ S 1 (t)- 


9Cin 

= 0 and m 4 = 0, this reads 
(A 2 + 12(7in (fW-fl 1 ^)) 

)*-*.* 

9(7in 


^A'| + I 2 A 222 (f 2 (t) — 0 2 (I)) 

) ! -a1 

9C*222 


(A 2 + 12A 222 (f 2 (t)-0 2 (t)) 

) 5 -4 

9C*222 


^A 2 + 12(7m (f x (£) — 0 4 (t)) 

)*-*# 


9(7111 


where 

(227) 


f(t) = f(xo-t). 

Inserting the last two lines into the first two, we obtain the following system for mj, to g, 
0 1 and 0 2 , 

mj = 2f* + 2 0 ! , 

(228) 3 _ . 9 _ * i ^ {1,2}. 

K? + 9C«i(f + 0 l ) = (A' 2 + 12 Ciuif - 0 1 )) 4 , 

In order to solve, for each i € {1, 2}, the second equation, we set 

Z t = (A 2 + 12Ciii(f — 0*)) ^, 

and notice that Z,, solves 

3 Zf + 4 Zf = 3(A 2 + UCiuff - 0 1 )) + 4 (AT? + 9 Cuiff + 0 *)) 

= 3 A 2 + 4Af + 72 C«iT, 

which is independent of 0 l . Denoting the right-hand side by d, i.e., 

(229) ci = 3 A 2 + 4Af + 72C„. 1 f, 

we see from the second equation in ( 12281 ) that 


( 230 ) 


Zjo - A? 

9 Cm 
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where Z lt o is the unique positive solution of 3 Z 4 + 4 Z 3 = a, which exists for any given 
Ci > 0 (i.e., for small enough |f |), since 3Z 4 + 4Z 3 is an increasing function of Z for positive 
Z. Notice that we have Zi t o = i/Kl when f = 0. Using u = u ,(y + U( r ) when x < 0, where 
U(i) and U( r ) are given by (1221D and (12241) . respectively, we have thus obtained an explicit 
solution in the vacuum region in terms of the initial profile f. Regarding then the crystal 
region, we obtain m x (x,t), i £ {1,2}, x > 0, from m l (x,t) = mo(to), where mj is given 
by the first line of (12281) . inserting the solution (1250)1 . and to is given implicitly, owing to 
(1217)1 . by the equation 

(231) x = m-o(to)^ 3 (1 - to)- 

This concludes the discussion of the decoupled case. 

5.7. A few words about the general case c ^ 0. We conclude this section with a 
few remarks on the generic case, where c ^ 0 and no Riemann invariants can be found. 
Obviously, the vacuum region can be treated in a similar way as in the decoupled case, 
the incident wave being represented by (122D) . and the ansatz (I223D for the reflected wave 
remaining valid, and leading to (12241) . Also, (12251) are to be used when applying the jump 
conditions (12201) . However, the crystal region cannot be dealt with as easily. Having no 
Riemann invariants at our disposal, we have to appeal to the evolution equations along 
the characteristics of the components w 1 = e* l -^u, i £ {1, 2, 3, 4}, of the spatial derivative 
of the solution u, as in the main body of the paper. Note from (fTOD) that 

Ai = y/m + R, A2 = Vm — R, A3 = —V m — R, A4 = —\Jm + R, 

so Ai, A2 are positive and bounded away from zero, while A3,A4 are negative and, again, 
bounded away from zero. The solution is trivial for x > Kit, the line x = Kit being a 
Ci characteristic in the crystal. A Ci or a C2 characteristic originating at a point in the 
non-trivial region x < Kit in the crystal intersects in the past the material boundary 
{x = 0}. Consequently, w 1 and w 2 must be given boundary conditions on {x = 0}. On 
the other hand, a C3 or a C4 characteristic originating at a point in the non-trivial region 
intersects in the past the line x = Kit which constitutes the boundary of the trivial region 
in the crystal. As a consequence, w 3 and w 4 satisfy trivial initial conditions oni = Kit. 
To deduce the boundary conditions for w 1 and w 2 on the material boundary, we impose 
the jump conditions (12251) . which, since the vectorfield is tangential to {x = 0}, imply 
similar conditions for the time derivatives: 


r d 


d 


d 


d 


=0, 


=0, 


=0, 

[dM 


In applying the first two of these, not only must we solve the system dnnp around the 
solution 

(Dy,D Z ) |( B#jB;t )_( 0]0 ) — (0>0)> 

but also its derivative with respect to t, namely the linear system 
f ^-Ey = Ki^- t D y + 6CinD y ^- t Dy + 2Ci 12 (^D z ^- t Dy + Dy^- t D?j +2Ci 22 D z ^- t D z , 

I |-E Z = K 2 ^-D z + 6C 22 2D z ^-D z + 2 C 122 (D y ^-D z + D z ^-D y ) + 2Cii 2 D y ^-D y . 
f at at at \ at at J at 

As in the decoupled case, the jump conditions not only provide the boundary conditions 

for w 1 and w 2 , but also give the boundary conditions for the two components (g 1 and 

g 2 ) of the reflected wave. Here, we use the fact that, by virtue of the equation (11811) . the 

components e* 1 of the time derivative of the solution u are given by: 

*id . i 

e —u = —XiW . 

at 

A key difference from the decoupled situation, aside from the characteristics inside the 
crystal not being necessarily straight lines, is that we expect the reflected wave to extend 
to all t > 0 at the crystal boundary (hence to all t > —x in the vacuum region). This is due 
to the fact that, inside the crystal, w 3 and w 4 will be non-zero throughout the non-trivial 
region x < Kit, generating non-zero boundary conditions at the material boundary for 
the reflected wave for all t > 0. This is because (see equations (ED) of Subsection 13.11) the 
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equations for and contain, respectively, the terms 27 3 2 w 1 w 2 and '2jf 2 w 1 w 2 

which are inhomogeneous from the perspective of the variables ( w 3 ,w 4 ). 
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